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Abstract 

New field theoretic renormalization group methods are developed to describe 
in a unified fashion the critical exponents of an m-fold Lifshitz point at the 
two-loop order in the anisotropic (m ^ d) and isotropic (m = d close to 8) sit- 
uations. The general theory is illustrated for the N -vector (ft 4 model describing 
a d-dimensional system. A new regularization and renormalization procedure 
is presented for both types of Lifshitz behavior. The anisotropic cases are for- 
mulated with two independent renormalization group transformations. The 
description of the isotropic behavior requires only one type of renormalization 
group transformation. We point out the conceptual advantages implicit in this 
picture and show how this framework is related to other previous renormal- 
ization group treatments for the Lifshitz problem. The Feynman diagrams of 
arbitrary loop-order can be performed analytically provided these integrals are 
considered to be homogeneous functions of the external momenta scales. The 
anisotropic universality class (N,d,m) reduces easily to the Ising-like (N,d) 
when m = 0. We show that the isotropic universality class (N,m) when m 
is close to 8 cannot be obtained from the anisotropic one in the limit d — » m 
near 8. The exponents for the uniaxial case d = 3, N = m = 1 are in good 
agreement with recent Monte Carlo simulations for the ANNNI model. 
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I. INTRODUCTION 



Formulated for the first time in 1975 by Hornreich, Luban and Shtrikman in the context 
of magnetic systems the Lifshitz critical behavior has encountered applications in many 
real physical systems. Some examples include high-T c superconductors ||f|,f|, ferroelectric 
liquid crystals || |7|, § , uniaxial ferroelectrics ||, some types of polymers [|l^,[ll],|12|,0] and 
magnetic materials Jl4, 15, 16, 17]. In particular, the confluence of a disordered, a uniformly 
ordered and a modulated ordered phase characterizes the special critical point associated to 
this critical behavior, known as the Lifshitz point. The modulated phase possesses a fixed 
equilibrium wavevector which vanishes continuously as the Lifshitz point is approached. 
When the components of this wavevector span an m-dimensional subspace, the system under 
consideration displays an m-fold Lifshitz critical behavior. When the order parameter has 
N components, and the space dimension of the system is d, the Lifshitz universality class 
is characterized by the set (N, d, m) . Whether m ^ d, the system presents the anisotropic 
Lifshitz critical behavior. Otherwise, the m = d case denotes the isotropic Lifshitz critical 
behavior. The isotropic case m = d near 8 can be treated, using similar theoretical tools, 
along the same lines of the anisotropic case. Thus, we shall focus our attention in these two 
types of critical behavior. 

In magnetic systems, the uniaxial (m = 1) Lifshitz behavior can be described by an 
axially next-nearest-neighbor Ising model (ANNNI) []18", 19|, which consists of a spin-1/2 



Ising model on a cubic lattice with nearest-neighbor interactions as well as next-nearest- 
neighbor antiferromagnetic couplings along one single lattice axis, the competing axis. The 
competition between the ferro and antiferromagnetic interactions in this system provokes 
a different critical behavior when compared to the pure Ising-like behavior. The magnetic 
compound MnP has been studied extensively in recent years, confirming the appearance of 



the uniaxial (m = 1) Lifshitz behavior which was obtained from theoretical |T^, [15], [16| as 
well as experimental |17| investigations. 



This model can be generalized by allowing the next-nearest-neighbor antiferromagnetic 
couplings along m directions, which represents a typical m-axial Lifshitz critical behavior. 
In case m ^ d, the system naturally admits two independent correlation lengths, namely ^2 
associated to spatial directions perpendicular to the competing axes, and £za associated to 
directions parallel to the m-dimensional competition subspace. At the Lifshitz point these 
two correlation lengths become related. In the isotropic behavior m = d close to 8, there is 
only one correlation length ^4. 

The field-theoretical representation of this model can be expressed in terms of a modified 
X(j) 4 theory containing higher derivative terms along the m-competing directions. It is given 
by the following bare Lagrangian density : 

L = \\ Vm 00 I 2 + \\ V(d-m) 00 T + 5 ~| Vm 00 ? + ^000 + A O0O- W 

The field theory treatment turns out to be simpler at the Lifshitz point, where T = Tl and 
5q = 0. In particular, the functional integral representation permits, at the Lifshitz point, 
the decoupling of the momentum integrals parallel and perpendicular to the competing 
axes. It would be interesting to find out whether this condition could make it possible the 
evaluation of Feynman diagrams to any desired order in a perturbative expansion. Then, the 
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critical properties of the system, like critical exponents, amplitude ratios and other universal 
amounts could be obtained analytically utilizing the renormalization group analysis along 
with e^-expansion methods. We shall consider this problem from a rather new perspective, 
which allows a solution in perturbation theory, and which provides an analytic tool that may 
prove useful in order to figure out the Lifshitz critical behavior in its complete generality 

In this work we present a detailed construction of the new renormalization group de- 
scription for the anisotropic and isotropic cases. This approach was inspired by an earlier 
suggestion made by Wilson JH] in order to obtain the critical exponents corresponding to 
correlations parallel or perpendicular to the competing axes in a manifestly independent 
manner. This framework was set forth in a previous letter |22|. We discuss the fundamental 
issues concerning this new renormalization group (RG) analysis leading to new scaling re- 
lations to the isotropic case, whereas in the anisotropic case it is shown that these relations 
are equivalent to previous scaling laws already derived. 

In the anisotropic behaviors, the existence of the correlation lengths ^ 2 and ^ 4 induces 
two independent characteristic external momenta scales which in turn are used to fix the 
renormalized field theory in the infrared regime. Since the theory is massless in the Lifshitz 
critical temperature Ti, the renormalized vertex parts have to be defined at nonvanishing 
external momenta scales. We denote by Ki the external momenta scale along the quadratic 
(noncompeting) (d — m) directions, whereas K2 is the external momenta scale along the 
quartic (competing) m-dimensional subspace. These external momenta scales originate two 
independent renormalization group flows in the parameter space. The renormalized coupling 
constants flow to two independent fixed points, depending whether the renormalization group 
transformation is over K\ or K2. At the loop order considered here they are shown to be 
the same, but it is suggested that this feature is preserved when the analysis is carried out 
for arbitrary loops. On the other hand, the isotropic case is characterized solely by the 
correlation length ^4. It induces only one characteristic external momenta scale, denoted 
here by K3 which is used to fix the renormalized vertex parts. 

Moreover, we calculate all the critical exponents at least at 0(e 2 L ) using a novel tech- 
nique of solving higher-loop Feynman diagrams inspired in this new renormalization group 
program. Our analysis is performed entirely in momentum space, which is particularly suit- 
able to tackle this problem. The Feynman diagrams are carried out with the help of a new 
approximation in the quartic momenta subspace, which is the most general approximation 
consistent with homogeneity. We will show how a former two-loop approximation presented 
earlier in the calculation of the critical exponents perpendicular to the competing axes [ PU| , |2B|] 
can be understood in terms of those calculated here using this more elaborate procedure. 
It is also shown that the relations among the correlation length exponents parallel and per- 
pendicular to the competing axes, namely via = \ v li, an d the anomalous dimensions of the 
fields, r)L2 = |??L4 are exact at the loop order considered in the present paper. This confirms 
the strong anisotropic scale invariance predicted before for this sort of system pl | . 

In Sec. II we set the formalism by defining the normalization conditions for the m-axial 
Lifshitz critical behavior for the anisotropic and isotropic criticalities. We show that two 
sets of normalization conditions can naturally describe the anisotropic situation without the 
need of introducing other dimensionful constants into the analysis. 

Section III contains the discussion of the new renormalization analysis for the anisotropic 
case in directions perpendicular to the competition clXGS, clS well as the new renormalization 
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group description along directions parallel to the competing subspace. 

The renormalization group treatment for the isotropic case is the subject of Sec. IV. The 
scaling laws are then obtained for this type of critical behavior. Since the scaling relations 
are different from those appearing in the anisotropic behavior, we point out that the isotropic 
and anisotropic behaviors are independent and cannot be obtained from each other. 

The evaluation of Feynman integrals is presented in Sec. V. We perform the one-, two-, 
and three-loop integrals using two different approximations. The first approximation in- 
troduced in |20],|23| is suitable to perform two- and three-loop integrals in order to obtain 
the critical exponents perpendicular to the competing axes, since it preserves the homo- 
geneity of the Feynman integrals in the external momenta perpendicular to the competing 
axes [ 25| . On the other hand, a new approximation is presented here which preserves the ho- 
mogeneity of the Feynman diagrams not only in the external momenta perpendicular to the 
competing axes, but also in the external momenta parallel to the competing m- dimensional 
subspace. Using a simple condition in the competing subspace, we calculated these integrals 
for arbitrary external momenta. 

In Sec. VI we calculate all the critical exponents for the anisotropic case using the scaling 
relations derived in Sec. III. It would be interesting to obtain the critical exponents using 
more than one renormalization condition in order to check their correctness. This is done in 
this section and in the following one. We also discuss our results comparing with alternative 
field theoretic treatments and new Monte Carlo simulations in d = 3 in the context of the 
ANNNI model (m = 1). 

Section VII presents the calculation of all the critical exponents for the isotropic case 
utilizing the new scaling relations obtained in Sec. IV. The analytical expressions obtained 
are new and to our knowledge are presented here for the first time. 

Finally, the conclusions are presented in Sec. VIII and further applications of the method 
described in this work are pointed out. 



II. NORMALIZATION CONDITIONS FOR THE LIFSHITZ CRITICAL 

BEHAVIOR 

From the bare Lagrangian given in (1) we can define renormalized quantities in terms of 
bare ones through the use of renormalization constants, or renormalization functions. Here 
we shall follow closely the standard A</> 4 field-theoretic approach. The interested reader 



should consult, for example, the Amit's book |26| or the original work by Brezin, Le Guillou 



and Zinn- Justin |27[| . These renormalization functions are fixed by the specification of the 

renormalization scheme used in order to define the renormalized theory. The renormalization 

functions are defined in terms of the renormalized reduced temperature and order parameter 

-i 

(magnetization in the context of magnetic systems) as t = Z^to, M = Z£ O and will 
depend on Feynman integrals. If the theory is renormalized at the critical temperature 
(t = 0), the infrared divergences instruct us to renormalize the theory in nonvanishing 
external momenta. Therefore, the renormalization constants at the critical temperatute Tl 
will depend on the external momenta scales involved in the renormalization program. 

We first consider the anisotropic behaviors. The Feynman integrals depend on two 
external momenta scales. We find convenient to define two sets of normalization conditions 
appropriate to calculate the critical exponents associated to correlations either perpendicular 
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to or along the competing axes p2 |. These external momenta scales were defined above to 
be Ki and k 2 , respectively. 

In order to make the calculation easier wherever more than one momentum remains finite, 
we choose the momenta at a symmetry point (SP). The normalization conditions which yield 
the critical exponents associated to correlations perpendicular to the competition axes are 
given by first setting all the external momenta along the competition axes to zero (k 2 = 0). 
Let pi be the external momenta perpendicular to the competition axes and associated to 
a generic 1PI vertex part. Then, the external momenta along the quadratic directions are 

2 

chosen as pi.pj = ^(4% — 1). This leads to (pi + Pj) 2 = k\ for i ^ j. The momentum scale 
of the two-point function is fixed through p 2 — k 2 — 1. Thus, we have the following set of 
renormalized 1PI vertex parts: 

rg ) (0,( 7l ) = 0, (2a) 
dY ( ${p,g 1 ). 

— 0^2 \p 2 =-i = l> ( 2b ) 

T { i\ P i,gi)\sp = gu (2c) 

(Pi,P2,P,gi)\sp = 1, (2d) 
rg' 2) (p, 9l )\ p 2 =Kl =0. (2e) 

Recall that the symmetry point implies that the insertion momentum in Eq.(2d) satisfies 
P 2 = (P1+P2) 2 = k\- 

The suitable normalization conditions to dealing with exponents along the competition 
axes are defined in a similar fashion. Firstly, one sets all the external momenta perpen- 
dicular to the competition axes to zero (k\ = 0). If k\ is the external momenta along the 
competition axes associated to a generic 1PI vertex part, the external momenta along the 
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quartic directions are chosen as k\.k'^ = ^f(4:5ij — 1). This implies that {k[ + k'^) 2 = k 2 for 
i ^ j. The momentum scale of the two-point function is fixed through k' 2 = k\ = 1. The 
analogous set of renormalized 1PI vertex parts is given by: 

rg ) (0,p 2 ) = 0, (3a) 
dk' 4 

r { *\k>,g 2 )\ S p = g 2 , (3c) 

{k'i,k 2 ,k',g 2 )\sp = 1, (3d) 

T$ 2 \k',g 2 )\ k/i=Kt =0. (3e) 

Note that, in principle, these two systems of normalization conditions seem to provide 
two renormalized coupling constants, which arise as a consequence of the two independent 
flow in the renormalization momenta scales K\ and k 2 . Apparently the whole description 
works with two coupling constants, namely g\ = iti(re?) 2 (and Ai = Moi(^i)~) associated 
to the flow in the momenta components perpendicular to the m-dimensional clXCS, clS well 
as g 2 = u 2 (k 2 )~ (and A 2 = ^02(^2)") associated to the flow in the momenta components 
parallel to the m-dimensional axes. Nevertheless, as will be shown, the situation simplifies 
at the fixed point: both couplings will flow to the same fixed point, at two- loop level, 



I K 4 



1, (3b) 
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indicating that this must be so in higher-loop calculations. The conceptual advantage is to 
treat independently the flow in the momenta along and perpendicular to the competition 
axes using these two coupling constants. Whether this can be done in a consistent manner 
is a separate problem, to be tackled in Sec. VI. 

The normalization conditions for the isotropic case (m — d near 8) can be defined 
analogously as those parallel to the competition axes for the anisotropic case. The symmetry 
point is chosen as follows. If k[ is the external momenta along the competition axes, the 

2 

external momenta along the quartic directions are chosen as k[.k'j = -j-(4^j — 1). This 
implies that {k[ + k'-) 2 = k\ for i ^ j. The momentum scale of the two-point function is 
fixed through k' 4 = k\ = 1. Then we have the following conditions: 



rf(0,s 3 ) = 0, (4a) 

Q^i \^=4 = ( 4b ) 

r ( £\K,g 3 )\ SP = g 3: (4c) 

T^\k[,k' 2 ,k',g 3 )\ s - P = l, (4d) 

rf 2 \k',g 3 )\ k , =4 =0. (4e) 



Notice that we have not mentioned the quadratic momenta scale K\ in the discussion of 
the isotropic behavior, for it is absent in this situation due to the Lifshitz condition 5 = 0. 

We can write all the renormalization functions and bare coupling constants in terms 
of the dimensionless couplings. Let the label r = 1,2,3 refer to the different external 
momenta scales involved in the general Lifshitz critical behavior, as discussed above for 
different normalization conditions in the anisotropic and isotropic cases. By expanding the 
dimensionless bare coupling constants u OT and the renormalization functions Z^ T ), Z^i^ = 
Z^t^Z^^ in terms of the dimensionless renormalized couplings u T up to two-loop order as 

u T = u T (l + a lT u T + a 2r i4)) (5a) 
Z^r) = 1 + b 2r ul + b 3r ul, (5b) 
Z02 (r ) = 1 + c lT u T + c 2T ul, (5c) 

along with dimensional regularization will be sufficient to determine all the critical expo- 
nents. 



III. RENORMALIZATION GROUP ANALYSIS FOR THE ANISOTROPIC CASE 

Given one bare theory, described by the Lagrangian (1), different versions of the renor- 
malized vertices can be constructed out of the original bare vertex parts. We shall explore 
now the freedom left in the definition of the renormalization momenta scales K\ and n 2 in 
the critical theory explained in the last section for the anisotropic case. 

We start by considering the renormalization group analysis along directions perpendic- 
ular to the competing axes. The renormalized theory is defined with only one quadratic 
nonvanishing external momenta scale K\. Let K x be the associated cutoff corresponding 
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to this subspace. The renormalized vertex parts for this case are defined in terms of the 
normalization constants and the bare vertices as : 



r 2(r?(P<(T). Q*(t), 9t, «r) = Z ^ T ) Z iHr)( r{N,L \Pi(r), Qi(r), K, K) (6) 

Sn,qSl,2^)\Q(t), Q(t), K, A t )|q2 ^2) 

where puj) (i — 1, N) are the external momenta associated to the vertex functions 
with N external legs and Qi( T ) (i = 1,...,L) are the external momenta associated to the 
L insertions of (f> 2 operators. We emphasize that pn\) (i = 1,...,N) refers to the external 
momenta components along the (d — m) dimensional subspace perpendicular to the com- 
petition axes, whereas p^2) are the external momenta components along the m-dimensional 
competing subspace. From our normalization conditions, it should be kept in mind that all 
quantities presenting a subscript r = 1(2) are calculated at zero external momenta compo- 
nents parallel (perpendicular) to the competing axes and are characterized by the momenta 
scale k\{k2). From the last section, uq t , Zm t \ and Zpw) are represented as power series 
in u T . The renormalization group invariance of the bare vertex with the momenta scale k t 
implies that : 

(^^)A.,A.[^7 ) ^ ) (r^ ) - s N , 8 L , 2 r\^)\ = o. (7) 

This in turn yields the following RG equations: 

Ll^ (T) (g T ,K T ))T { ^(p i(T) ,Q i{T) ,\ T ,K T ) = 5 N)0 5 Li2 (k; 2t )^B t , (8) 

where B T is a constant used to renormalize and 

- dg T 

Pr(gr, Kr) = (^t^ 1 )a t ,A t (9a) 
7#r)(gT,«r) = {K>T dln Q* {T) )\ T ,K T (9b) 

dlnZ^^ 

1<P{t)K9t,Kt) = («T q k )\r,Ar- (9c) 

are functions of g T and k t only, though they are functions of A r implicitly. Notice that 
r^' 2 ' ) is different from all other vertices since the RG equation presents an inhomogeneous 
term in the left hand side due to its additive renormalization. We shall treat this additively 
renormalized vertex part later on. The above expressions correspond to the limit A r — > oo, 
which are naturally finite, even if A T , ^( T ) and Z^( r \ diverge at this limit. It is worth 
expressing all of these quantities in terms of dimensionless bare and renormalized coupling 
constants. We now turn our attention to discuss the central issue of the new dimensional 
considerations which will be useful for the subsequent dimensional analysis. 

Consider the volume element in momentum space for calculating an arbitrary Feyn- 
man integral. It is given by d d ~ m qd m k, where q represents a (d — m)-dimensional vector 
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perpendicular to the competing axes and k denotes an m-dimensional vector along the 
competing subspace, respectively. The Lifshitz condition 5q = suppresses the quadratic 
part of the momentum along the competition axes. Nevertheless, there is still a contri- 
bution from the quartic momenta contained in the inverse critical (t = 0) free propagator 
1 (q, k) = (k 2 ) 2 + q 2 . In order to be dimensionally consistent, the canonical dimension 
in mass units of both terms in the propagator should be equal. There are two ways out of 
this outstanding situation. 

The former idea, inspired in Ref. is to introduce a dimensionful constant a in front 
of the first term in the Lagrangian (1), along with its renormalization function Z a , as was 



done by Mergulhao and Carneiro ||28|| . This idea implies that the momenta scales parallel or 
perpendicular to the competition directions play the same role in this discussion and there 
is only one coupling constant. Denoting the components of the quartic external momenta 
with a subscript a and the quadratic components with a subscript (3, they set the following 
renormalization conditions: 

dT R (K -k,<r,g,K) 



dk 2 lk ? =K2 



d^jk, -k,a,g,K) 
dki 

,(4), 



7fci=K 2 



r ( R > (k i ,a,g,K)\ SPa = g, 

(h,k2,P,<T,g,K)\sp a = !i 
r S' 2) (p, -P, v, 9, K ) |op*=« a = 0, 

where the SP a means a 2 was chosen at zero quadratic external mo- 

menta. This choice of renormalization points makes the renormalization constants a inde- 
pendent as claimed by those authors [p9j] . However, a is still a relevant length (momentum) 



scale and this fact should be reflected on its dependence in some normalization constants. 
Therefore, starting with a dimensionful a parameter and making it dimensionless in the 
end of the calculation as they chose does not seem to be consistent, since the quartic and 
quadratic momenta scale play the same role and have the same canonical dimensions in 
this approach. Notice that the last 4 of these equations together with the critical theory 
condition on the renormalized two-point vertex part naturally defines an independent set of 
normalization conditions along the competing axes. In fact, if the quartic momenta is rede- 
fined through ah* = k£ such that a is absorbed in the new quartic momenta, this implies 
that a^ki a kj a = k' ia kj a = k' 2 ( 4< ^~ 1 ) with k' 7^ k. Then, one has 5 normalization conditions 
along the quartic subspace as described in the last section. 

On the other hand, the first of these equations is calculated at ak\ = 0. Intuitively it 
should be complemented with 4 more normalization conditions with nonvanishing external 
quadratic momenta perpendicular to the competing subspace. This is what was done in 
the last section for directions perpendicular to the competing axes. Thus, if we have two 
different momenta scales k and k' and setting them equal is equivalent to have Mergulhao 
and Carneiro's renormalization conditions, with 5 more normalization conditions along the 
quadratic directions. Thus if one trades a by an independent external quartic momenta 
scale k', it still remains the 5 extra normalization conditions which in their approach are 
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undefined. Nevertheless, they recovered the former anisotropic scaling relations |T[] using 
this reasoning. They used their symmetry point in order to treat the cases m = 2, 6 in the 
context of an e^-expansion |[29|| . 

There is an alternative based in a recently proposed method which does not use the di- 
mensional constant a but allows the realization of a dimensional redefinition of the momenta 
components along the quartic competing subspace |22| . This later view inspires the subse- 



quent discussion and shall be used throughout this paper. Let [q\ = M be the mass dimension 
of the quadratic momenta. The consistency of the Lagrangian density (1) on dimensional 
grounds requires that [k] = Ma. This is equivalent to perform a dimensional redefinition of 
the momenta along the competing long as the condition 5o = is satisfied. The vol- 

ume element in momentum space d d ~ m qd m k has mass dimension [d d - m qd m k] = M d ~^ . The 
dimension of the field is obtained by requiring that the volume integral of the Lagrangian 
density (1) is dimensionless in mass units. It follows that [(f)] = M^ d ~^~ l . 

The iV-point Green function can be expressed dimensionally as [G^ N \xi, ...,xn)] = 
[(f)] N = M~( d ~^~ N . The associated one particle irreducible (1PI) vertex parts have dimen- 
sion in mass units [r^^)] = [G^ N \x t )}[V]- N [G (2 \x t )}~ N = M^ d ~^ +N . In momentum 
space, the Fourier transform is obtained by integrating over each one of the coordinates. 

Removing the momentum conserving 5-function, we have [r^(fcj)] = M N+<yd ~^ 2 . 

As usual, the exponent of M in the above relations is called the canonical dimension of 
the quantity. If the physical quantity O has canonical dimension [O] = M A , then under a 
transformation of the length scale associated to M — > olM it implies that O = a A 0. There- 
fore, all dimensionfull parameters are transformed under a transformation in the lengths 
(or external momenta). Hence, it is useful to describe the theory in terms of dimensionless 
parameters. As the coupling constants are associated to r^, we can write g T = u t (k 2t ) ^ , 
and A r = u 0t (k 2 J)~^ , where e L = 4 + y — d. 

In terms of the dimensionless couplings defined above, the renormalization group equa- 
tion can be rewritten as: 

d d 1 e 

^ T d/T + ^ T ~du ~ 2^ 7<A(r) ^ + L ^ 2 w( Mr ^ r S? = S N,oh,2(K-; 2r )^B T (u T ) } (11) 

and from now on we can forget about the cutoff's A T , bearing in mind, however, that they 
should be kept fixed in all stages of the analysis. The functions 

Pt = (Kt 1 T : ), (12a) 



dlnZ^ T ) 



l^M = (3 T ~—^ (12b) 
dlnZ^ (r) 

7<A 2 (r)(«r) = ~Pr ^ (12C) 

are calculated at fixed bare coupling A T . The /3 T -functions can be cast in a more useful form 
in terms of dimensionless quantities, namely, 

Pr = -re L (-^-) • (13) 
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Note that the beta function corresponding to the flow in K2 has a factor of 2 compared to 
that associated to the flow in K\. As usual, they are power series in u T , with coefficients 
which depend on e L . Let us analyse the simplest case L = 0. The solution can be expressed 
in terms of characteristics. The characteristic equation is given by two independent flows in 
the coupling constants induced by the flows in the momenta scale K\ and K2, i-e, 

p T ^^ = P(u T (p T )), (14) 
dpr 

with the initial condition u T (p T = 1) = u T . Using the characteristic equation for u T we can 
change variables from a variable x T to u T , through the relation: 

f pT ft t \\ dx T f uApT) fM 

J I X T Ju T Pt{U t ) 

Thus, small values of x T in the left-hand-side correspond to the neighborhood of the zero's 
of (3 T in the right-hand-side. For the anisotropic case, the solution to the renormalization 
group equation reflects the two-parameters group of invariance, and can be expressed in the 
form 

T{ R(l)( k i( T ),u T ,K T ) = exp[-^- f ^ {T) (u T (p T ))—]T { ^ ( l ) (k i{T) ,u T (p T ),K T p T ). (16) 

From our dimensional analysis, the dimensional redefinition of the momenta along the 
competing axes results in an effective space dimension for the anisotropic case, i.e., (d— y). 

Thus, we find the following behavior for the 1PI vertex parts under flows in the external 
momenta: 

r (N) , , , r(iV+(d-f)-^X2) r N fPr dX T 
r R(T)(Prh(T),U T ,K T ) = p T exp[-— 7«(t)(Ut(Pt)) J (17) 

^iZ(T)(^»( T )' U t(Pt), K T p T ). 

It is helpful to present the explicit expressions for vertex parts calculated either at zero 
quartic external momenta or at vanishing quadratic external momenta. The renormalized 
vertex parts calculated at zero quartic external momenta is given by 



r R{ [ ) (p l k i{1) ,u 1 ,K 1 ) = Pl exp[-— 7^(1) ] (18) 

r S^])(^(i) 5 «i(pi)>«ipi)- 



The dependence of the renormalized vertex parts is quadratic in the external momenta 
perpendicular to the competing axes. Therefore, the analysis is completely similar to the 
pure A</> 4 -theory, with the replacement — > e. From this analysis, we can identify the label 
t = 1 with the subscript L2. Then, we could have written 7^^) = 70(^2) and 7^2(1) = 702 (z,2)- 
On the other hand, the renormalized vertex parts at zero quadratic external momenta 
can be expressed as 

(N) 2(N+(d- — )-^—Zl) N j'P' 2 dX2 

^R(2)(p2ki(2),u 2 ,K 2 ) = p 2 2 2 exp[-— ^ 70(2)(M 2 (p 2 ))— ] (19) 

T R{2)( k i(2)^ U 2{p2)^2p2)- 
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The difference is that the canonical dimension is twice as big as the canonical dimension of 
the vertex parts calculated at zero quartic momenta. Then, we can make the identifications 
7^(2) = 7</>(L4) and 7^2 ( 2 ) = 7^,2 ( L4 ). The last equations imply that a change in the external 
momenta scale is equivalent to the multiplication of the vertex function by that scale to the 
power of the canonical dimension of the function, followed by a modified coupling constant, 
which flows with the characteristic equation, and an additional factor. 

It is interesting to analyse the vertex functions at the infrared fixed points, since this will 
determine the scaling laws and the critical exponents associated to correlations perpendicular 
and parallel to the m-dimensional competing subspace. The analysis can be carried out 
by assuming that there are two independent fixed points, defined by (3 T (u*) = 0. The 
renormalization group equation leads to a simple scaling property at the fixed points. It 
implies the following solution to the vertex functions: 

C(M, W ,«.> T ) = ^*>-^*HFW«> (20) 

For N = 2, we have 

rg T) (p T fc (T) , <, Kr) = ^'XfeX. «r)- (21) 

The quantity ^^)(u* T ) can be interpreted in the following way. If the field theory is free, a 
change in the external momenta scale will produce a change in the free vertex which 
scales with the canonical dimension of the vertex, that is 

i(r) {Prkiir)) - pr 1 (t) (fcj( T ) J . (22) 

We then define the dimension of the field (ft as 

rSy(p.^ )) = P ; [(d -- ) -^ 1 r^(A; l ), (23) 

such that in the free theory <i° = — j 2 — 1 is the naive dimension of the field. At the fixed 
point, the naive dimension is modified due to the presence of interactions, such that the 

cZ— — 

nontrivial effect is the appearance of the anomalous dimension, i.e., d^ T ) = — j 2 — 1 + 
When N = 2, this identifies the anomalous dimensions of the field in the anisotropic 
situation, namely that associated to the change in the external momenta scale perpendicular 
to the competing axes 771 = r]L2 = 7<Ki)( M *) as we ^ as ^ ne other corresponding to the change 
in the external momenta parallel to the competition subspace = VL4 = 7</>(2)(' u 2)- 

This can be easily generalized to include L insertions of (f> 2 operators in quite an analogous 
way, such that the RG equations at the fixed point lead to the solution ((TV, L) ^ (0, 2)) : 

(N,L) 1 , * \ _ T[ly + (d ~-> 2 lL \ 2 — +jL V(r) 



V r(t) (P k i(r) , PPi(r) , <, Kt) = P 

-i(AT,L 
R(r) 

Thus, if we write at the fixed point 



X ffl ) ft(T).ft(T). ,( >r)' ( 24 ) 



11 



T-i(iV,L) / , * x r[(d-f )-Nd t f,}+Ld !j> 2 (N,L) n * x , _\ 

T R(r) yPh{T),PPi{r),U T ,K T ) = p T T' r{t) (k i{r) , p i(r) , U T , K T ) , (25) 

the anomalous dimensions of the insertions of 2 operators are given by = —2t+^^2^{u* t ) 
and will be related to the critical exponents v\ = 1/12 and V2 = ^L4 as we shall see in a 
moment. 

In order to find the scaling relations we must go away from the Lifshitz critical temper- 
ature (t 7^ 0) staying, however, at the critical region 5 = J35|. Above the Lifshitz critical 
temperature, the renormalized vertex parts for t 7^ can be expanded as a power series 
in t around those renormalized vertices at t — 0, provided N 7^ 0. We can now apply the 
differential operators 

d d 1 d 

to r^-j (ki(r)) such that we find 

00 t L d d 1 

Orr^Cfyr),*, <,«,-) = E^|[ K -^r + ^^-- 2 iV7 ^ (r) ^ T ' ) + L702(r) ^ ^ 

The result is that each term in the sum vanishes because of the RGE for 

r i?(r^(^(r))Pj(r), «*, «r)- Hence, we obtain 



T-r^r + ^ T ^7" ~ T^^mC^t) + V(r)( M r)^]r^ } (A; i(r) ,t,<,K T ) = 0. (28) 



The solution is a homogeneous function of the product of fcj( r ) (to some power) and £ solely 
at the fixed point u*. As the value of u T is fixed at u*, we shall omit it from the notation of 
this section from now on. It is given by: 



r^l } {k i[r) ,t,K T ) = K7^F^\k l(T) ,K T f^y). (29) 

If we define 9 T = — 7^2 ( T ), and using dimensional analysis, we find 

1 i?(r) I'M"*") > ' ""rj — A'r /Vr 

FgPOvV).^" 1 ^)^*)^). (30) 

By choosing p T = v(-2t) 9t+2t , and replacing back in (29), the vertex function depends 
only on the combination k^ T ^ T apart from a power of 1 Since the correlation lengths £ T are 
proportional to t~ UT , we can identify the critical exponents v T as 

^ 1 = 2r + ^ = 2r- 7 ; 2(T) . (31) 

According to our conventions, these equations are equivalent to the following scaling relations 
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^ = 2-7^(1), (32a) 

^ = 4-7^(2)- (32b) 
As a matter of convenience, we could have defined alternatively the function 

x , („ \ - R dln ( z <P{r) z <t>{T)) hql 

1^{t){Ut) - -Pr 7^ • {ii) 

In that case we would have found the equivalent formulae 

v ll = 2~Vl2- 7>(i)K), (34a) 

v ll = 4 - Via - %2 {2) . (34b) 

Hence, at the fixed point all correlation functions (not including composite operators) scale 

at T > Tl, since they are functions of ki( T )^T only. For iV = 2 we choose p T = kr T ), 
the external momenta. Then r^ T ^(fc( T ), t, k t ) = k 2r ~ rh k1 t /(fc( T )^ r ). The critical situation 
is characterized when £ r — > oo and fc( T ) — > such that /(fc( r )£ T ) — > Constant. Defining 
/r = (^(r)?r) 2r_, ' T /(A:( T )^ T ), we have 

rw^t,^) = (fc (T) e r ) 2r -^^/r(fc( T) ? r ). (35) 

(2) —1 

The susceptibility is proportional to t~ lT as fc( T ) — > 0. Thus, since = X( r )> we can 
identify the susceptibility critical exponents 

TV = i/ T (2r - ?7 T ). (36) 

These relations are equivalent to the relations: 

7^2 = ^L2(2 - 77 L2 ), (37) 



7i4 = ^L4(4 - r/ L4 ). (38) 



(0,2) 



The specific heat exponents can be obtained by analysing the RG equation for T^ T ^ 
above Tl at the fixed point, i.e. 



;«r^: + 7^ (T )(2 + ^)) r K = (0*^(0, (39) 

where B T (u*) is given by 

(^)^S r «) = -Z| 2(T) K T A r [;f(Q (T) ;-Q (r) ,A r )|^ )=K? . (40) 

It is a inhomogeneous part which has no dependence in the reduced temperature t. The 
bare vertex function r?'-^ is calculated as before in the limit A T — > oo, with a fixed bare 
coupling constant, which renders B T (u* T ) finite in this limit when (d — y) = 4. This renor- 
malized vertex part consists of the addition of the homogeneous (temperature dependent) 
and inhomogeneous pieces. The general discussion given so far for the vertex part r^y^ 
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will be useful to determine the homogeneous part of the solution. Indeed, at the fixed point 
the obvious generalization of the solution for is given by: 



Therefore, the temperature dependent homogeneous part for h will scale at the fixed 
point as: 



r { Ril[ h (Q(r),-Q(r),t,K T ) = (Q (r) , -Q (t) , K T t "^M ) . (42) 

This will be identified with the specific heat at zero external momentum insertion Q( r ) = 0. 
Using the results of our dimensional analysis 

r$;i fc )(Q(T),-Q ( T),*,/0 = ^^^^^^^(^Qm.-p^^w.p^p; 1 ^), (43) 



and replacing this into the solution at the fixed point, we find 

rg?>*«M, -Q M ,t, Kr) = (44) 



-1 



x F( ' 2 )(p r - 1 Q (r) ,-p r - 1 Q (T) ,p; 1 /c T (p r - 2 ^)V M ). 

Again we choose p T = k t (4f) * t+2t • Substituting this in last equation, taking the limit 
<5( T ) — > and identifying the power of £ with the specific heat exponent a T , we find: 

Tfl 

a T = 2-r{d--)v T . (45) 

Let us analyse the inhomogeneous part. First, take <3( r ) = 0. Second, choose a particular 
solution of the form: 

C p (u T ) = {k 2 t t )^C p {u t ). (46) 
Replace this into the RG equation for at the fixed point. Then, it is easy to obtain 

w) = ( ^ rT(d : f) _ 2 ^(<). (47) 

Collecting both terms we have the following general solution at the fixed point: 

= (v-")*(C^)- + VTT{d \ } _ 2 BA<)). (48) 

Let us describe the situation for T < T^. It can be illustrated for the case of magnetic 
systems. The renormalized equation of state relates the renormalized (1PI one-point vertex 
part) magnetic field with the renormalized vertex parts for t < via a power series in the 
magnetization M. One has: 
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OO -y 

H (T) (t, M, U T , Kr) = £ JJi MNr R{r)( k i(r) = ^ *, U r, «r), (49) 

where the zero momentum limit must be taken after realizing the summation. The magnetic 
field satisfies the following RG equation: 

(*tq£ + Ptq£ ~ 2 N 7«r)( N + M dM ] + ^ t m )H ^ M ' w - ^) = °' ( 5 °) 
At the fixed point we have the following form for the equation of state: 



H( T) (t,M,Kl) = Kr 2 /il T (fi T M*,M V(T) )' ( 51 ) 

Once again, we use dimensional analysis arguments to obtain the following expression under 
a flow in the momenta : 

H (T) (t,M,K T )=p T T [ ^ +1] H (T) ( — , - *L). (52) 

Pt M— ^ — i] 

We choose p T to be a power of M such that: 

Pr = Kr[ f[(rf M f) - 2] ] T|(d ^ 2) - 21+ - ; (53) 

and from the scaling form of the equation of state Hr T )(t, M) = M &T f( — ^-), we obtain the 
remaining scaling laws after the identifications Si = Sl2, Pi = Pl2, S2 = Sla, P2 = Pla, '■ 

(d - + 2 - r, L2 
S " = ((( _ ?) _ 2 + , u . < 54 »> 

! 777 

P L 2^^L2((d--)-2 + r, L2 ), (54b) 

2(^-^+4-^4 
= 2(d-f)-4 + W (54C) 

I 777 

/?L4 = 2^( 2 ( d " 7) " 4 + %4), (54d) 

(54e) 

which imply the Widom ^ L2 = Pl2{Sl2 — 1) an d Rushbrook a i2 + 2/?£2 + 7l2 = 2 relations 
for directions perpendicular to the competition axes. These relations are also valid for 
directions along the competing axes. So far, the effect of considering this new dimensional 
role played by the momenta scale along the competing quartic subspace, together with the 
definitions of the critical theories either at vanishing quartic or quadratic external momenta 
have induced two independent set of scaling relations for the critical exponents. Nevertheless, 
when performing the diagramatic perturbative expansion, we shall find out that some of 
these exponents are not independent. 

There is one curious fact relating these scaling relations and those first obtained by 
Hornreich et al. |]. K v la = \ V L2 and 772,4 = 2r/£ 2 to all orders in perturbation theory, the 
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In 


= 1L2 


= 1L, 




= Pl2 


= Pl, 




= 5 L 2 


= Sl, 


a L4 


= a L2 


= 0L L 



hyperscaling (Josephson) relations are the same in either formulation. In the formulation 
presented here this feature will be found in Sec. VI from the perturbative analysis up to the 
loop order considered in this paper. In addition, one obtains 

(55a) 
(55b) 
(55c) 
(55d) 

In that case, there is a complete equivalence among the scaling relations in both formulations. 
It is worth to emphasize that the advantage of the approach presented in this section is 
the splitting of the scaling laws into independent renormalization group flows parallel and 
perpendicular to the competition axes. Then, instead of claiming that the two momenta 
scales corresponding to the components perpendicular and parallel to the competing axes 
are equal 0, the most important conclusion in our approach with two coupling constants is 
that they will flow to the same fixed point in the critical regime as will be shown later. 

We can make a comparison among our scaling relations below Tl with the ones obtained 
by Mergulhao and Carneiro. In their work, they defined the space dimension of the Lifshitz 
system D = 4f + dp, where d a = m and dp = d — m (see Eq.(25) in fl28| ). Therefore, 
D = (d— y)> the same effective space dimension as ours. 

The trouble is in the introduction of a. The normalization conditions defined through 



Eqs.(20)-(24) of Ref. [28[] mix a with the two external momenta scales in a nontrivial way. 
They recognized, however, that the normalization conditions they defined are independent 
of a. Thus if one makes the choice 7* = 0, Eq. (58) in Ref. is just the same as the one 
obtained here for Sl, and Eq.(60) in Ref. |28| is equal to that obtained here for (3 L . The 
Eqs. (46)-(51) in their article taken together with 7* = yields trivially vl± = \ V L2 and 
77^4 = 2r] L2 . Last but not least, if one takes the bare a dimensionless, as was done by those 
authors, the whole argument is invalidated since its dimensionfulness was assured from the 
beginning of the discussion in the regulation of the free critical propagator. Introducing 
a is consistent provided it is considered as a dimensionful ammount in all stages of the 
calculation. In other words, a is not required at all, since the flow in a can be absorbed in 
the quartic momenta scale using our dimensional redefinition. 



IV. RENORMALIZATION GROUP FOR THE ISOTROPIC BEHAVIOR 

The procedure to analyse the isotropic case is quite analogous. Some care must be taken. 
Whenever r appears as a subscript, like in a quantity A T , one sets r = 3 in order to be 
consistent with the notation employed so far. The dimensional analysis is a bit different. 
The volume element in momentum space is again d d ~ m qd m k. Whenever d = m, the volume 
element is now d m k. As before, [k] = M2. Accordingly, the volume element has dimension 
[d m k] = M~ . The dimension of the field in mass units is [(f)] = M~~ x . When the conserving 
(^-function is removed the 1PI vertex parts have dimensions [r^(fcj)] = M n+! %~ n! t . Then, 
make the continuation m = 8— e^. The coupling constant has dimension A3 = = M~ . 

In terms of dimensionless quantities, one has the renormalized = ^3(^3)^ and bare 
^3 = ^03(^3)^ coupling constants, respectively. Again, the functions 
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Ps = («3^) (56a) 
7«3)(«s)=A^^ (56b) 
7*(s)(«s) = -/5a^^ (56c) 

are calculated, as before, at fixed bare coupling A3. The beta function can be expressed 
in terms of dimensionless quantities as f3$ = — £l{ 91 q^ One should notice that the 
beta function for the isotropic case does not possess the overall factor of 2 present in the 
anisotropic beta function obtained from renormalization group transformations over the 
quartic momenta scale K2- This feature is a strong nonperturbative suggestion that the 
isotropic critical observables cannot be obtained from the anisotropic ones and vice-versa. 

Under a flow in the quartic momenta, from our dimensional analysis, the dimensional 
redefinition of the momenta along the competing axes results in an effective space dimension 
for the isotropic case, i.e., (y). Thus, we find the following behavior for the 1PI vertex parts 
r^(3) under a flow in the external momenta: 

r£g)(P3^3,«3) = pl [N+f ~ Nf] exp[-f JF 70(3)Mp 3 ))ff] (57) 

X ri R(l)(h,U 3 {p 3 ),K 3 p 3 ). 

Notice that we put aside the notation k^ T ), etc., used in the anisotropic analysis in favor 
of hi, etc. since there is only one quartic momenta scale in the isotropic case. At the fixed 
point, we also have a simple scaling property for the vertex parts r^x, namely: 

r%te(fc,,u| lK 3) = ^ +? -" ¥, -^"» w) (58) 

For N = 2, we have 

rg; 3) (M, «S, K3) = P3 _7 ' (3)K) r§3)(A;, «$, « s ). (59) 

We can now identify i] L/L = t] 3 = 7^(3) (ifjj) as the anomalous dimension for the isotropic 
case. This is the analogue of the analysis we performed for the anisotropic case. In the free 
theory — -| — 1 is the naive dimension of the field. At the isotropic fixed point, the naive 

m 

dimension is modified due to the presence of interactions, such that d^ — -| — 1 + The 
generalization to include L insertions of (f> 2 operators is quite straightforward and can be 
written at the fixed point as ((N, L) 7^ (0, 2)) : 

T R(3)(p3 k i,p3Pi,U 3 ,K 3 )=p 3 T^^PiiU^Ks). (60) 

Thus, if we write at the fixed point 

T R(3){p3h,P3Pi,U 3 ,K 3 )=p 3 r K m '{ki,p i: U 3: K 3 ) : (61) 
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the anomalous dimension of the insertions of <fi 2 operators is given by = —4 + 7^(3) (uQ. 

Above the Lifshitz critical temperature, the renormalized vertex parts for t 7^ can be 
expanded as a power series in t around those renormalized vertices at t — 0, provided N ^ 0. 
We can now apply the differential operators 

d d 1 d 

3 = k 3 — + fc— - -N lm (u 3 ) + 7 ^ (3) ( U3 )t— (62) 

to r^g^fcj, t, U3, /C3). The mechanism is similar to that discussed in the anisotropic cases. 



This vertex part is a power series on t, with each individual coefficient vanishing by making 
use of the RGE for T ^^\ki,pi, u%, k 3 ). Then, we find 



d d 1 d 

^ 3 l^ + l33 du 3 ~ ~ 2 Nlm M +7^(3)(«3)*^]r2j3 ) (fci,t,u^K 3 ) = 0. (63) 

The solution is a homogeneous function of the product of fcj (to some power) and i only at 
the fixed point u%. The solution reads: 

Nl 4>(3) -1 

rgg ) (fc i(S ),t,«S,«s) = K^F^iki^t^m). (64) 

All the exponents generated by the renormalization group flow along the scale n 3 will be 
denoted by a corresponding L4 subscript. If we define 9 3 = — 7^2(3), one can use dimensional 
analysis to obtain 

r (N) /, , x _ 2[N+f -f f ]-f VL4 f VL4 
1 R(3)l Ki ' r ' K 3j — P3 K 3 

x F^\p^k t ,(p^K 3 )(p^). (65) 

One can choose p 3 = k 3 {\) e »+ 4 , and replacing it in (64), the vertex function depends 
only on the combination k£ L4 apart from a power of t As ^4 is proportional to t~ ULA the 
critical exponent v 3 = v L4 can be identified as 

^4 1 = 4 + ^ = 4-7; 2(3) . (66) 
Again it is convenient to define the function 

7<^ 2 (3)(%) = -As 7^ — — ■ (67) 

Then, one can easily find the following relation 

v ll =4-Vl4~ 7^(3) (u$) ■ (68) 

At the fixed point all correlation functions (not including composite operators) scale at T > 
T L) since they are functions of k^ L4 only. For = 2 we choose p 3 = k, the external momenta. 
The two-point vertex part can be written in the form r^L (k, t, k 3 ) = k 4 ~ VL4 kI la f '(k^Li) ■ The 
main point is that when ^4 — ■> 00 and k — > 0, simultaneously, then f{k^Li) —> Constant. 
Defining f 3 = (k^ L4 ) 4 -^f(k^ L4 ), we have T^ 3) (k, t, « 3 ) = (A&4) 4 -"" ^3" / 3 (^l 4 ) • The 
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susceptibility is proportional to t 7i4 as fcj — > 0. Since = x 1 , the susceptibility critical 
exponent is given by 

7L4 = ^L4(4 - r] L4 ). (69) 

The scaling relation for the specific heat exponent can be found from the RG equation 
for r^gj above T L at the fixed point, namely 

(*s^ + 7fc(8)(2 + t^mm = (^B 3 (ul), (70) 

where 

(n^B 3 (u* 3 ) = -4(3) K 3^r^(Q; -Q, As), (71) 

is the inhomogeneous part which does not depend on t. Recall that the bare vertex function 
r^ 2 - 1 is calculated as before in the limit A 3 — > oo, with a fixed bare coupling constant, which 
renders -83(143) finite in this limit when m — 8. This renormalized vertex part is made out of 
the addition of the homogeneous (temperature dependent) and inhomogeneous pieces. The 
general discussion given so far for the vertex part is helpful to obtain the homogeneous 

part of the solution. At the fixed point we have the following generalization of the solution 
for T (N ' L) - 

R(3) 



r^ffe, Q i} t, k 3 ) = 4 ^4 N > L \ Pl , Q u M Vo) ). (72) 



The temperature dependent homogeneous part for V^B. h scales at the fixed point, i.e., 

r { R§, h (Q,-Q,t,m) = ^ 2i;2(3) F^ 2 \Q,-Q,K 3 t^). (73) 

This vertex function is to be identified with the specific heat at zero external momentum 
insertion Q = 0. Using our dimensional analysis one finds: 

rggj >fc (Q, -Q, t, « s ) = P^+^vr^tfQ, -tfQ, PsX ^ 1 « s ). (74) 

Substitution of this equation into the solution at the fixed point yields 

, , /i4 2[f-4]+2 7 * 2 ^ -2 7 * 

R(3),h ! 



rS 3 2 L(Q,-Q,^ 3 ) = P 3 2 W W ^^(^-^Q.ft 1 ^^) 7 ^)). (75) 



We make the choice p 3 = ft3(4r) e 3+ 4 . Replacing this into the last equation, taking the limit 
Q — > and identifying the power of i with the specific heat exponent «L4, we find: 

a LA = 2 - mz/ L4 . (76) 

The description of the inhomogeneous part is as follows. First, take Q = 0. Then, choose a 
particular solution of the form: 
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C p (u 3 ) = (4)^C p (u 3 ). (77) 
Now replace this into the RG equation for T^gj at the fixed point. Therefore, one gets to 



C„(«S) = (^-^^B^)). (78) 
The general solution at the fixed point is just the sum of the two pieces, and is given by 

rgsl = (^(^r" + TT^SW))- (79) 

We now turn our attention to analyse the scaling relations when the system is below 
the Lifshitz critical temperature T < Tl. The renormalized magnetic field is related to the 
renormalized vertex parts for t < and the magnetization M through 

OO 1 

H {3) (t,M,u 3 ,K 3 )= 2-^M Jv r^(fc i = 0;t,«3,«3), (80) 

N=l ly ■ 

where the zero momentum limit must be taken after realizing the summation. The magnetic 
field satisfies the RG equation: 

(«3^ + - -N 7m (u 3 )(N + M m ) + l4>H3) t-)H {3) (t, M, ui, «i) = 0. (81) 

The equation of state at the fixed point reads: 



H {3) (t, M, k 3 ) = k 3 2 h 3 (K 3 M>iL4 , «!tV(3) ). (82) 

Dimensional analysis arguments lead to the following expression under a flow in the external 
momenta: 

H {3) (t,M,K 3 ) = pp +1] # 3 (4 ^ttt,-). (83) 

P3 p 3 U 1 P3 

The flow parameter p 3 is chosen to be a power of M such that: 

M 

P3 = ^ 3 [ nf ^]] m - 4+ "«, (84) 

and from the scaling form of the equation of state H^(t,M) = M &LA f\ — V- ), we obtain 
the following scaling laws: 

m + 4 - r] L4 

Ola = 7- , (85a) 

m - 4 + r] L4 

Pla = (m - 4 + rj L4 ) , (85b) 
which imply the Widom 7^4 = I3l 4 (5l 4 — 1) and Rushbrook a^ 4 + 2f3 L4 + 7^4 = 2 relations. 
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The scaling relations for the anisotropic case Eqs.(lla),(12) in Ref. [JJ for d = m are 
consistent with the isotropic case. Note, however, that this cannot be given a rigorous 
meaning, for the appearance of v L2 and r] L2 in the equality in Eq.(llb) of Ref. [|IJ invalidates 
the argument for the isotropic case as these exponents are no longer meaningful. Notice that 
the impossibility of finding scaling relations for the isotropic case in the original framework 
is due to the lack of the independent flow in the external momenta scale K3 along the 
quartic subspace. In the early treatment JTJ, the quartic momenta was not independent 
to be varied freely, but was fixed from the variation of the quadratic scale. Without its 
free variation, which is possible since this quartic term in the propagator does not have 
the same canonical dimension as the quadratic one, no renormalization group flows along 
the competing directions can be defined whatsoever. Thus, this new renormalization group 
method permits to go further in determining the Lifshitz critical universal properties of the 
system for arbitrary m. 



V. THE EVALUATION OF FEYNMAN INTEGRALS 

In order to calculate universal quantities like critical exponents, we must calculate some 
Feynman integrals. The perturbative loop expansion shall be our starting point with the 
ei = 4 + y — d being the perturbation parameter for the anisotropic situation. For the 
isotropic case, the perturbation parameter is = 8 — m. 

We have to express the solution of the Feynman diagrams in terms of e^, resulting in 
the e^-expansion for the universal critical ammounts. Again, there is also a very important 
difference among the anisotropic and isotropic behaviors. From a technical viewpoint, the 
anisotropic behaviors present two types of integration along the two momenta subspaces, 
whereas in the isotropic situation there is only one subspace to be integrated over. We shall 
treat them separately. 

The anisotropic behavior is described using two different approximations. We shall briefly 
discuss the first analytical approximation developed for evaluating higher-order Feynman 
diagrams which are needed in the calculation of the critical exponents perpendicular to the 
competing axes for the anisotropic Lifshitz behavior. It points out the necessity of some 
sort of condition among the quartic loop momenta in different subdiagrams, leading to the 
homogeneity of the integrals in the quadratic external momenta scales. We employ the set 
of normalization conditions with vanishing quartic external momenta as described in Sec. II. 
This piece of work was done in collaboration with L.C. de Albuquerque and the details can 
be found in [p0| ,|23|. 



Nevertheless, with the renormalization group description presented here, this approxima- 
tion is not sufficient to describe the critical exponents along the competing axes. It does not 
yield the solution of the integrals as a homogeneous function of both quadratic and quartic 
external momenta scales yet. The former approximation described above is then generalized 
to obtain the solution of the integrals for arbitrary quadratic and quartic external momenta 
scales. Using the new interpretation for the momenta scale along the quartic direction given 
in the last two sections, the calculation of these integrals is not a complicated task, provided 
a certain condition among the quartic momenta is fulfilled. With this new technique all the 
critical exponents in the anisotropic cases are obtained. This picture can be considered the 
main result of the present work. 
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The isotropic behavior can be developed along the same lines of the latter approach to 
the anisotropic case. The condition among the quartic momenta is also required in order 
to guarantee homogeneity of the Feynman integrals in the quartic external momenta scale. 
The new approximation is sufficient to complete the unified analytical description of the 
Lifshitz critical behavior in its full generality, at least at the loop order considered here as 
will be shown in this section. 

In order to verify the renormalization scheme independence of the critical exponents, it 
would be interesting to obtain the critical exponents using more than one renormalization 
procedure. In fact, as will be proven later, the use of normalization conditions or minimal 
subtraction of dimensional poles yield the same critical exponents. Thus, we shall present 
the results in the most appropriate form for calculating the critical exponents in these two 
renormalization prescriptions. 

A. Anisotropic 

In order to calculate universal quantities like critical exponents, we must calculate some 
Feynman integrals. We start by listing all the relevant integrals which are necessary to find 
out the critical exponents. These integrals are 



/ 



d d - m qd m k 



i((k + K'yy + (q + P) 2 ]((k 2 ) 2 + qi) 
where I 2 is the one-loop integral contributing to the four-point function, 

d d - m q 1 d d - m q 2 d m k 1 d m k 2 



(86) 



r (r-" l qi<r-" l q 2 d" l k l d" l k 2 

J S? + (^i) 2 ) (5 + m 2 ) [fa +Q2+ py + iih + h+ k') 2 ) 2 } ' (8,) 



is the two-loop "sunset" Feynman diagram of the two-point function, 

d d ~ m q x d d - m q 2 d m k x d m k 2 



n = 



{ql + {klY) ((p - qi y + ((K> - hf ) 2 ) {ql + {klY) 

x ? . (88) 

{qi-q2+P3) 2 + {{ki-k 2 + k' 3 Yf 

is one of the two- loop graphs which will contribute to the fixed-point, and 

d d - m q l d d - m q 2 d d - m q z d m k 1 d m k 2 d m k z 



I 



(q 2 i + (^i) 2 ) (Ql + (^l) 2 ) (Ql + (^l) 2 ) [(0i + 02 - P) 2 + {{h + k 2 - k'Y) z ] 

x 5 (89) 

(0i + 03-p) 2 + ((£:i + A;3-A;') 2 ) 2 



is the three-loop diagram contributing to the two-point vertex function. Now we proceed 
to calculate these integrals using two different approximation schemes. The philosophy to 
be addopted is to simplify the calculation by making use of the homogeneity hypothesis, as 
shall become clear in the following subsections. 
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1. The "dissipative" approximation 



As this approximation is only suited to calculate the integral as a function of the quadratic 
external momenta, we set the external momenta at the quartic directions equal to zero, i.e. 
k' = k[ = k' 2 = k' 3 = 0, and K' = k[ + k' 2 . We shall use dimensional regularization for the 
calculation of the Feynman diagrams. 

Let us find out the one- loop integral I 2 . With our choice of the symmetry point, and 
introducing two Schwinger's parameters we obtain for I 2 : 



d d - m qd m k 



((P) 2 + (3 + P) 2 )((P) 2 + g 2 ) 

x J d d ~ m q exp(-(a a + a 2 )q 2 - 2a 2 q.P - a 2 P 2 ). 



roo roo / r 

J J rfairfa 2 (y d m k exp(— (ai + a 2 )(£; 2 ) 2 } 



(90) 



The q integral can be performed to give 



J d d - m q exp(-(oti + a 2 )q 2 - 2a 2 q.P - a 2 P 2 ) 

1 ^,d—m., d-m OL\OL 2 P 2 > 

= -S d - m T( — —){a 1 + a 2 ) 2 exp( ■ 

z 2 cti + a 2 



(91) 



For the k integral we perform the change of variables r 2 = k\ + ... + k^. Now take z = r 4 



The integral turns out to be: 

J d m k exp(-( ai + a 2 )(A: 2 ) 2 ) = (ls m )T{j){ ai + a 2 ) 
Using Eqs. (90) and (91), I 2 reads 



(92) 



2 - "*V4 

r°° f°° OL\d 2 P 2 
x / / da.\da 2 exp( ; ) [a\ + a 2 ) 



o Jo 



d rn_ 

.2 4 / 



(93) 



The remaining parametric integrals can be solved by a change of variables followed by a 
rescaling The integral is proportional to (P 2 ) 2 . Now we can set P 2 = k 2 = 1. Using 
the identity 



r(a + fox) = r(a) 1 + b x ip(a) + 0(x 2 ) 



(94) 



where if)(z) = -^lnT(z), one is able to perform the e^-expansion when the Gamma functions 
have non integer arguments. Altogether, the final result for I 2 is: 



-s m s d ^ m T(2 - -^)r(— ) 



1 



1 + [*a]m £l 



(95) 



where [i 2 ] m = 1 + |(^(1) — ^{2 — r j)). The factor inside the brackets in Eq. (94) is absorbed 
in a redefinition of the coupling constant. Then the redefined integral is: 
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/2 = ^-(l + N m ei) • (96) 

Note that this expression involves no approximation. This simple result is a conse- 
quence of the absence of the zero quartic external momenta. Had we considered it from the 
beggining, we would have obtained an intermediate integral that could not be integrated 
analytically. We shall discuss this issue later in the next subsections. 

However, when we go on to calculate higher loop integrals, some sort of approximation 
is required, since these integrals are complicated by the fact that even with zero external 
quartic momenta, the quartic loop momenta mix up in different subdiagrams in a extremely 
non-trivial form. As an example, we discuss 1%. It is given by: 

j (p K>) = f d d - m q 1 d d - m q 2 d m k 1 d m k 2 



Setting K' = 0, the integral can be evaluated as outlined in PUJ. Before making our 
approximation, one can choose to integrate first either over the loop momenta (qi,ki) or 
over (q 2 , k 2 ). The loop integrals to be integrated first are referred to as the internal bubbles. 
By solving the integral over q 2 first, we obtain 

1 v( d-m f d d - m qi d m h 

I" f°° d ai da 2 ( ai + a 2 ) zi ^exp{ ^-(qi + pf) I cT^e^ 1 ^ ! e -^((fci+* 2 ) 2 ) 2 ! . ( 98 ) 

Jo Jo ot\ + a 2 J 

Now we can consider the approximation. In order to integrate over k 2 , we have to 
expand the argument of the last exponential. This will produce a complicated function of 
ai,a 2 ,ki and k 2 which cannot be integrated analytically. Considering the remaining terms 
as a damping factor to the integrand, the maximum of the integrand will be either at k\ = 
or at k\ = —2k 2 . The most general choice k\ = —ak 2 yields a hypergeometric function. 
The choice k\ = —2k 2 implies that k\ varies in the internal bubble, but not arbitrarily. 
Its variation, however, is dominated by k 2 through this constraint, which eliminates the 
dependence on k\ in the internal bubble. At these values of k\, the integration over k 2 
produces a simple factor to the parametric integral proportional to (a% + a 2 )~~ . This 
allows one to perform the remaining parametric integrals in a simple way. After performing 
these integrals, they produce the factor ((gi + P) ) a - . Note that the diagrams J 3 and 
I 5 contributing to the two-point function receive the factor tt-Ku- after integrating over the 

Z 4 

quadratic momenta in the external bubble. This factor will not be present in the isotropic 
case, since there is no integration over quadratic momenta to be done in this case. The 
resulting solution to Is(P, 0) is a homogeneous function of the external momenta P, given 
by: 



8 — m €l 



(99) 



The implementation of this constraint on higher-loop integrals proceeds analogously. The 
constraint turns all these integrals into homogeneous functions of the external quadratic 
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momenta scale. One can then choose the symmetry point as P — 
order to define the renormalized vertices via normalization conditions. 



1, for example, in 



Using this constraint we can easily find the following results at the symmetry point |20 

2 4(l + 3[ 



The integrals I' 3 and I5 are given by 



8 — m €l 



1 + [i 2 ]m + 



4 _ m. 

^ 2 



(100) 



(101) 



3(2 - f ) 4 



1 + 2 [i 2 ] m + 



cy m 

Z 4 



(102) 



Note that the leading singularities for J 2 , 1± are the same as their analogous integrals 
in the pure 4 theory. However, I' 3 and I' 5 do not have the same leading singularities for 
they include a factor of ^3™) • We t nen introduce a weight factor for Jg and namely 

(1 — y), so that they have the same leading singularities as in the pure (f) A theory. The 
main drawback of this approximation is the failure to treat the isotropic case. Furthermore, 
the introduction of weight factors to the two and three-loop diagrams is rather undesirable. 
Moreover, the constraint among the quartic loop momenta does not allow loop momentum 
conservation along the quartic subspace in higher-loop diagrams. It is then appropriate to 
name this approximation "the dissipative approximation". 

It is obvious that some important detail is missing. A proper solution of the Feynman 
integrals should be expressed as a homogeneous function of both external momenta scales. 
We proceed to discuss the novel approximation which presents this property. 



2. The orthogonal approximation 

Before considering the integrals to be performed, let us derive some useful formulas which 
relate Gamma functions with certain intermediate parametric integrals. They will allow us 
to define a new analytic dimensional regularization procedure in the competing subspace. 

The simple integral 

-il 1 

expi^—ax^dx = a » — IY— ), (103) 
/i 11 

can be generalized to the m-sphere. We shall analyse the case // = In. Take r 2 = x\ + ...+x^ n . 
After that take z = r 2n . Thus, 

/oo I /"OO r 

dx\...dx m ex r p{—a[x\ + ... + x'^ n ) n ) = — / dzexp(— az)z^~ 1 / dVt m . (104) 
-00 ' 2n Jo J 

The angular integral will produce the area of the m-dimensional sphere, yielding 

dx x ...dx m exp(-{xl + ... + x 2 m ) n ) = ^-T(^)a^S m . (105) 

In In 
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One can write this identity in a different way. After choosing r 2 = x\ + ... + x 2 m take y = r 2 
and as the integral is given by the expression above, we obtain the intermediate result: 

/ dyy-~ x exp{-ay n ) = -a~ T{— ). (106) 
J o n 2n 

Henceafter we shall keep n — 2. The following step is to calculate the integral 

/OO /'OO 
exp( — ax A — bx 2 )dx = 2 / exp(— ax 4 — bx 2 )dx. (107) 
-oo JO 

The exact answer is given in terms of a Bessel function of a certain combination of a and b. 
We wish to pick out only the piece which yields the correct homogenous function of a, i.e., 
only one term of the series. This can be achieved as follows. Choose y = x 2 . One obtains 

/oo ^2 „oo ^ _ 1 
exp(— ax 4 — bx 2 )dx = exp( — ) / exp(—a(y H ) 2 )y~dy. (108) 
-oo 4a Jo 2a 

We then choose y' = y + implying that 

r oo exp(-ax 4 - te 2 )d:r = exp[^] [f °° exp(-ay' 2 )(y' - ±)^dy> 

b 



^exp{-ay> 2 ){y' -D^dy'}. (109) 

Since we are dealing with convergent integrals, we can perform the approximation (y' — 
7^)~ = y'^z + and the remaining terms will be subtracted from the last integral, which is 
a sort of error function. The original integral is then approximated by its leading contribution 

/oo £)2 /-oo 1 b 2 \ i 1 

exp(—ax 4 — bx 2 )dx = exp(—) / exp(—ay' 2 )y'^^dy' = exp(—)-T(-)a~ 1 . (HO) 
-oo \n Jo 4a 2 4 

It can be shown in a straightforward way that for the m-sphere this result generalizes to 

f-oo exp[-a(x\ + ... + x 2 m f - b{x\ + ... + x 2 m )\dx x ...dx m exp(^)S m x 

^ exp(-ay' 2 )y^- 4 dy> = exp{^)\S m Y^ )a"f . (Ill) 

We now focus our attention in the integral 



i 



oo \{x\ + ... + x 2 m ) 2 + 2a(xf + ... + x 2 m ) + m 



21/3 • 



;ii2) 



Take r 2 = x\ + ... + x 2 m . Make the change of variables in the radial coordinate z = r 2 . After 
that take z' = z + a. We then obtain 

f {(xl+...+xl l )2+2a(x 2 1 +...+x2 n )+rn 2 }0 = \ Sm 

m _ i m i 

r roo (z — a)"2" <iz ra (z — a) "2" rfz i /-■ -■ q\ 

[JO (z'2 +m '2)/3 JO (z'2 +m '2)/3 J) l 11 ^ 

where m' 2 = m 2 + a 2 . Taking z" = z' 2 , expanding the numerator in the first integral, i.e., 
keeping only the leading term and getting rid of the infinite terms to be subtracted from 
the second integral, one can write this integral in the approximated form 
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J [(x-(+...+x^y 2 +2a(xf+...+x'i l )+m^ ~ 4°™l m " > 

r <?sp>. (iM) 



We have all ingredients to perform Feynman integrals for arbitrary m. We start by consid- 
ering the simplest integral, the one- loop integral contributing to the coupling constant, that 
is, 

I 2 = / ./ '■■"""■ k . (115) 

J [((k + K') 2 ) 2 + (q + P) 2 ]((k 2 ) 2 + q 2 ) V ; 

We can use two Schwinger parameters and integrate over the quadratic momenta. Using 
the formula derived above 



/ 



exp(-p 2 )A = ^r(^), (116) 



we obtain 



1 r( d-m y°° [°° aia 2 P\ 
h = ^d-nJ-{ — r — )/ / da 1 da 2 exp( ■ ) 

Z Z JO JO Oil + <*2 

x( ai + a 2 )-(^) |«exp(-a 1 (A; 2 ) 2 -a 2 ((A; + K') 2 ) 2 ). (117) 

We can now expand the argument of the last exponential. This integral cannot be performed 
analytically. We are interested in the solution of this integral in a form that preserves 
homogeneity in both external momenta. Some simplifying condition should be tried to 
achieve this goal. 

The most general approximation to calculate this type of integral which is homogeneous 
in the external momenta scales can be understood as follows. In the first place, if we set 
k.K = inside the integral, that has the virtue of eliminating odd powers of the quartic 
external momenta. Thus, the integral becomes: 

J d m kexp(- ai {k 2 ) 2 - a 2 {{k + K') 2 ) 2 ) = J d m k 
exp(-( ai + a 2 )(A: 2 ) 2 - 2a 2 k 2 (K') 2 - a 2 ((K') 2 ) 2 ). (118) 

Using Eq.(lll), we have for the last quartic momenta integral: 

/ «exp(-a 1 (A; 2 ) 2 - a 2 ((k + K') 2 ) 2 ) = S m \T{^ + a 2 )~T 

x exp(-^f^). (119) 

We can then express the integral in the following form: 

r r**** exp(- a ' a ' (P2 + ((g ' )2)2 ) («, + <*)-(!-*). (120) 

JO JO Oi\ + o> 2 



X 



Take x = a^P 2 + (K' 2 ) 2 ) and y = a 2 (P 2 + (K 12 ) 2 ). After that, define v = Thus, the 
parametric integrals can be done easily by this change of variables. Then, use the identity 
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T(a + bx) =T(a) 1 + bx ip{a) + 0{x 2 ) , (121) 



where ip(z) = 4-lnT(z). This will result in the following expression for I 2 : 



1 ri 

2 



[\s {d . m) s m v{2 - f )r(f )](i - f v(2 - f ))r(f ) 

x J, 1 dv(v(l - v)(P 2 + ((K') 2 ) 2 ))^ ■ (122) 

This is a homogeneous function (with the same homogeneity degree) in (P, K') just as 
advertised. But this is not the answer yet. The factor [^S^- m )S^r(2 — ^)T(^)] can be 
absorbed in a redefinition of the coupling constant. Hence, we shall absorb exactly this 
factor after performing each loop integral. Furthermore, the last integral can be expanded 

as 

£ dv(v(l - v)(P 2 + {{K'ff))^ = 1 - jL(P, K% (123) 

where 

UP, K') = C dv ln[v(l - v)(P 2 + ((K') 2 ) 2 )]. (124) 
Jo 

Thus, we find the following result for this integral: 

l2(P,K') = j-{l + ([ l2 ] m -l)e L - € fL(P,K'^ . (125) 

This is the form suitable for renormalizing using minimal subtraction. On the other hand, 
for normalization conditions one has: 

/2Sft=/2Sft = ^(l + Nm£L) , (126) 

since L(SP l = SP 2 ) = -2, with SP 1 = (P 2 = 1,K' = 0) and SP 2 = (P = 0,(K') 2 = 
1). When we calculated I 2 (P, K' = 0) in the last subsection, the orthogonality condition 
k.K' = between the loop momenta and the external momenta along the quartic subspace 
was trivial. In the calculation of I 2 (P, K'), the orthogonality condition allowed the solution 
to this integral with the correct homogeneous properties in both external momenta scales. 

We can now turn our attention to the higher-loop integrals. The simplifying condition 
k.K' = for the one-loop integral can be easily generalized to the higher-loop integrals 
by stating that the loop momenta in a given bubble (subdiagram) is orthogonal to all loop 
momenta not belonging to that bubble. Let us see how this works in the calculation of the 
"sunset" two-loop integral I 3 contributing to the two-point function, given by the following 
expression: 



j (p R n = f d d ~ m q x d d ' m q 2 d m k x d m k 2 

31 ' ] J {ql + {k 2 ) 2 ) (<1 2 2 + (k 2 2 ) 2 ) [fe + q 2 + P) 2 + ((fc + k 2 + K'Yf] ' 



We can choose to integrate first over the loop momenta (q± , k\) or over (q 2 , k 2 ). The loop 
integrals to be integrated first are referred to as the internal bubbles. By solving the integral 
over q 2 first, we obtain 
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1 d — m.rd m q 1 d m k 1 r°° r°° -(d- m ) 

I 3 (P,K) = -S d _ m r(—)J ql + {kl)2 J q J q da ld a 2 (a 1 + a 2 ) * 

( qi + P) 2 ) [ d m k 2 e-^ ik ' )2 e^ k ^ +K ' )2)2 . (128) 
«i + a 2 J 



OL\OL 2 

exp 



We define K" = k\ + K' into the argument of the last exponential, and integrate over k 2 
using the condition k 2 .K" = 0. Make the change of variables k\ = p and integrate over k 2 
(or p). Using Eq.(lll) we find: 

I" r daida 2 (ai + a 2 )-^ + ^exp{-^^[{ qi + Pf + {{k x + K') 2 ) 2 ]). (129) 
The parametric integrals can be solved as before and we have: 

We can now use Eq.(94) and absorbing the angular geometric factor for the first loop integral 
we obtain: 

UP, K') — — (1 + [i 2 ] m t L ) / « r . (131) 

Let is(P, K') be the last integral above. In the remainder, we employ a Feynman parameter. 
The integral is(P, K') can be expressed in the following form: 

^K') = ^±^J o 1 dxx^ 



r d d - m qi d m h 

J (a? 4 2.rP ^ 4- tP^ 4- M - rMhlV 4 .rfYfc, 4- /TWU+^ ' 



After that, take the orthogonality condition fci-K 7 = 0. In order to solve the integral over 
the quadratic momenta q\ we shall make use of the relation: 

f d d - gi i r(^)r(a-^)(m 2 -^-% 
7 (gf + 2fc. 9l + m 2 )« ~ 2 r(a) d ' m ' { ' 

Thus, we obtain 

■ / P ^ ^- m r(2-f + C£ )r(-i + g + 6L ) ,1 ^ 

«s(p,^) = ^ ^ y o ^ 

/ ^ (134) 

7 ((k 2 y + P^x(l-x)+x(2k 2 K' 2 + (K' 2 ) 2 ))- 1+: f + ^' y ] 

Now, using Eq.(114), we can integrate over the quartic momenta k\ obtaining: 
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s d - m s m r(2 - f + f )r(f )r(-i + e L ) 



8 T(f) 

x C dxx^-^xil -x)(P 2 + (K 12 ) 2 )} 1 -^. (135) 
Jo 

Expanding the resulting T functions and absorbing the geometric angular factor discussed 
above, we find: 

h(P, K') = (P 2 + (K' 2 ) 2 )^{\ + 2[i 2 ] m t L - -e L - 2e L L 3 (P, K')\ (136) 



where 



L 3 (P, K') = f 1 dx(l - x)ln[{P 2 + {K' 2 f)x{l - x)\. (137) 
Jo 



At the symmetry points SP\ or SP 2 , it can be rewritten as 



hsp, = hsp 2 = — ^(1 + 2{i 2 } m e L + -e L ). (138) 

o€l 4 

From the above equation we can derive the expressions: 

W= ^) = 4p 2 (- ^) = 8^(1 + 2[< 2 ] m CL + ^). (139) 

Let us now proceed to discuss the other required loop integrals. Consider 

d d - m q 1 d d - m q 2 d d ~ m q 3 d m k 1 d m k 2 d m k 3 



I 



(ql + (k!) 2 ) {ql + (k 2 2 ) 2 ) (q! + (kj) 2 ) [( Ql + q 2 - P) 2 + ((h + k 2 - K'YY] 

x (140) 

(gi + g 3 - P? + m + k 3 - K'Yf 1 ; 



which is the three-loop diagram contributing to the two-point vertex function. This integral 
is symmetric under a change in the dummy loop momenta q 2 — > q 3 and k 2 — > k 3 . Let us 
analyse the integrations over q 2 , k 2 and q 3 , k 3 . We use the condition k 2 .(k\ — K') = when 
integrating over k 2 as well as k 3 .(ki — K') = when performing the integral over k 3 . The two 
internal bubbles, which are represented by the integrals over (q 2 , k 2 ) and (q 3 , k 3 ), respectively, 
give actually the same result, namely I 2 ((qi — P), (ki — K 1 )). Next take P — > — P, K' — > —K'. 
Therefore, we obtain the following expression: 

i „ r jd—m n f ]mt, 

up, m = ^ + 2 M . tI ) / M - M - - ^- - g , )am . (hi) 

We employ a Feynman parameter in analogy to what was done in the calculation of I 3 and 
working out the details we find: 



h(P, K>) = (P 2 + (K' 2 f)-l(l + 3[ i2 ] m e L -e L - 3e L L 3 (P, K')), (142) 



At the symmetry points SP±, SP 2 we find: 
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hs Pl {= -Qpr) = hsp 2 (= -q^jt) = g^"(l + 3[i 2 ] m e L + -e L ). (143) 



We are left with the task of calculating one of the two-loop diagrams contributing to the 
four-point function 

d d - m q 1 d d - m q 2 d m k 1 d m k 2 



(ql + (k 2 ,) 2 ) {(P - Qi y + {{K> - h) 2 ) J (q 2 + (k 2 ) 2 ) 

1 



(<?i -Q2+ P3) 2 + ((fci -k 2 + k : 



(144) 



Notice that P = pi + p 2 , Pi (i = 1, •••,3) are external momenta perpendicular to the com- 
peting axes whereas ZT' = + fc 2 , and k\ (i = 1, 3) are the external momenta along the 
competition directions. We can integrate first over the bubble (q 2 ,k 2 ). Using Schwinger 
parameters, and absorbing the geometric angular factor for the first bubble we obtain: 



1 /■ d d - m q x d m k x 

J4 ' ' 



;i + M m eL)/ ((? 2 + (A ,2 )2)(( p 



l-J^J (g2 + (A a )a) (( p _ gi)2 + _ fci)2)2) 

x " ■ (145) 

[(91 + ps) 2 + m + 

Using a Feynman parameter one can write this in the form 

"i t d d - m q x d m k x 



fm^Lj J dz J 



[ql - 2zP. qi + zP 2 + (kf) 2 + z((K' 2 ) 2 + 2K' 2 k 2 )} 



2 



x ir, (146) 

[(Qi + Ps) 2 + ((h + k',) 2 ) 2 }^ 



where we defined the quantity / m (e^) = ^-(1 + [i 2 ] m ^L)i which is the one-loop subdiagram 
with the angular factor already absorbed. Using another Feynman parameter to fold the 
two denominators in the last expression, integrating over pi, k\ (recalling the orthogonality 
condition already stated), the integral turns out to be 



\Ue L ) MEM Si d yy(l- y)¥^ J* dz 



yz(l-yz)(P 2 + (K' 2 ) 2 ) 



+ 2/(1 - y)(pi + mm + 2yz(l - y)( P3 .P + (k' 3 ) 2 (K>) 2 ) 



Sd-mSm- (147) 



The integral over y is singular at y — 1 when = 0. We only need to replace the value 
y = 1 inside the integral over z [P6fl , and integrate over y afterwards, obtaining after the 
absorption of the geometric factor 



c 2 



h = ^(l + 2 [i 2 ] m e L - -e L - e L L(P, K')). (148) 



This is the most appropriate form to carry out the renormalization using minimal sub- 
traction. In terms of normalization conditions, we find the value of this integral at the 
symmetry points discussed before: 
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(149) 



Thus, we have successfully devised a new regularization procedure to calculate Feynman 
integrals whose propagators have quartic powers of momenta. It is tempting to define the 
measure of the m-dimensional sphere in terms of a half integer measure. In fact, taking 
k = p 2n , one has d m k = d^p = -^p^~ l dpdVL m . Hence, the approximation ammounts to 
take the new "measure" d^p invariant under translations p' = p + a. 

Note that this approximation is much better than the former dissipative approximation 
for the following reasons. First, we do not have to introduce any diagram factor, since 
after absorbing the geometric angular factor the leading singularities in are the same as 
those in the standard (j) 4 field theory present in e (m — 0). Second, we have an expression 
in terms of arbitrary external momenta, which permits the computation of all the critical 
exponents in a completely independent manner using renormalization group transformations 
either perpendicular or parallel to the competition axes. Third, this can be easily adapted 
to the isotropic behavior. The three weak points of the dissipative approximation have been 
overcome in the orthogonal approximation. 



B. Isotropic 

The new orthogonal approximation can now be used to obtain the solutions to the 
Feynman integrals in the isotropic case. At the Lifshitz point 5q = all the quadratic 
momenta disappear and only quartic momenta are defined. 

Consider the integral 

r d m k , s 

h = / o • (150) 

J ((k + K ) 2 ) (k 2 ) 2 V ; 

It is the isotropic counterpart of the one-loop integral contributing to the four-point function. 
Using two Schwinger parameters and the orthogonality condition k.K' = 0, we find: 

/f OO f OO 
d m k j j da^e-^^^^- 2 ^'^^-^^ 2 . (151) 

Now performing the transformation k 2 = p the volume element transforms to d m k = 
^p^^dpdQm = d~p. The former quartic integral turns into a quadratic integral over 

m 

p. After neglecting the infinite terms which change the measure d~k under the translation 
y' — y + ^, only the leading contribution is sorted out and we have 

J ne - a{ ef-^ = J d f pe -a P *- bp s a -f e £ hQ)S m . (152) 

Replacing this result into the expression of I 2 {K'), one finds 

I 2 (K') = ]r(^)S m r r da 1 da 2 (a 1 + a 2 )"* exp(-^^). (153) 
4 4 Jo Jo a\ + a 2 

Now, use a change of variables and a rescaling to realize the remaining parametric integrals 
analogously to what was done in the anisotropic case. Making the continuation m = 8 — e L , 
the integral can be expressed in the following intermediate step: 
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I 2 (K') = -T(2-^)TC-f)S m J o dv[v(l-v)(K a )*\-'+, (154) 
and the integration over v produces the result 

/ 2 (K') = — (1 + L(K' 2 ))]. (155) 

We absorb the factor of S m in this integral through a redefinition of the coupling constant. 
Henceafter we shall absorb this factor when calculating each loop integral in analogy to the 
discussion for the anisotropic behavior. Note that this absorption factor is different from the 
one appearing for the anisotropic case in the limit d — > m = 8 — ex- In the anisotropic case 
the geometric angular factor becomes singular in the above isotropic limit indicating the 
failure of this attempt of extrapolating from one case to another. This is a more compelling 
technical reason for the statement that the isotropic and anisotropic cases are completely 
distinct. Therefore, 

I 2 (K') = ±-[l-%(l + L(K a ))]. (156) 

In terms of a symmetry point (K' 2 ) 2 = 1 convenient whenever normalization conditions are 
used, we obtain 

h(K') = -[I + (157) 

We can go on to evaluate the other required higher loop integrals. The systematics is the 
same: solve the subdiagrams using the intermediate step Eq.(153) and then use Feynman 
parameters to solve the parametric integrals left over. 

Let us calculate 



J ((k 1 + k 2 + K') 2 ) 2 (k 2 ) 2 (k 2 ) 2 ' 



the two-loop "sunset" Feynman diagram of the two-point function in the isotropic case. 
Integrate first over k 2 . Take K" — ki + K' and use the condition k 2 .K" = to obtain: 

^ + fi/ K( , 1 + ^ w - (159) 

Now using a Feynman parameter, integrating over k±, taking m = 8 — ex, and employing 
the formula 



/ 



d^q „ 1 m )r(a - f )(m 2 - k 2 )f~ a S m 



(q 2 + 2k.q + m 2 ) a 4 T(a) 
one can express I3 as 



(160) 



gmF( VnV + ° ) [1 + T ] f Q dx[(K' 2 ) 2 X (l- X )}^^. (161) 
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We can rewrite this expression in terms of the integral L 3 (K' 2 ) defined in the last section in 
the form 

_ S m T(2-f)T(-l + f) e L 1 3e L e L 
h ~ 4^f(f) [1 + T ]l 2 " To 7 " ~2 L3{K )] " (162) 

Expanding the Gamma functions and absorbing S m , it is easy to show that 

h = - [ - T ^-[l + e L (--L 3 (K' 2 ))]. (163) 
At the symmetry point, this can be expressed as: 

1 + (164) 



16e L L 8 
leading to 

We can carry out the calculation of the other integrals using the same reasoning. The 3-loop 
integral J 5 is given by 

d m hd m k 2 d m k 3 (166) 



((*! + k 2 + iT) 2 ) 2 ^ + fc 3 + iT) 2 ) 2 (A; 2 ) 2 (fc 2 ) 2 (A; 2 ) 2 ' 

where we took for convenience K' — > —if'. The integrals over A; 2 and k 3 are identical. Hence 

^5 = ^-[l + y]/— — r .,, 2 iii n2 „ • ( 167 ) 
e L 2 J [((fci + if) ) ] 2 (A; 2 ) 2 



Employing the same techniques as in the calculation of I 3 we obtain 

h = -\t^^ + ^{-~ -MK' 2 ))]. (168) 



At the symmetry point, the derivative of I§ with respect to the external momenta is given 
by 

QT 1 

g^lm = ii=-^[l+«d. (169) 

The two-loop graph contributing to the 4-point function in the isotropic situation is 
j = [ d™hd m k 2 

1 J (kimx* - h) 2 ) 2 (k 2 2 ) ah - k 2 + k> 3 ) 2 f ' 

where K' = k[ + k' 2 . It can be integrated using the orthogonal approximation following the 
same steps of the calculation of the anisotropic counterpart. We simply quote the result 
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h(K") = ^[l~^(l + 2L(K"))]. (171) 

L 



At the symmetry point it is given by: 



It is worth to point out that the integrals 1% and 1$ have not the same leading singularities 
as in the usual <fi 4 . Therefore, any attempt to use the counterterms of the usual <fi 4 theory 
would lead to erroneous results for the critical exponents in the isotropic case. 

A recent calculation of the critical exponents 77^4 and ula for the isotropic case was 
presented by Diehl and Shpot who fixed all the leading singularities equal to those appearing 
in the standard (f> 4 theory ||30|| . Moreover they "choose" the following angular factor: 

1 , „r(5- |)r 2 (| -2) , s 

~ ^4 lira) 

If one sets d = 8 — 6l, whenever the Feynman integral under consideration presents a double 
pole in €l, this angular factor will give contribution to the simple poles in e^. This happens 
for example in the two-loop integrals if and I4. Then, their calculation of these critical 
exponents cannot be trusted, since further terms were not taken into account in the 
evaluation of those integrals. It seems that this factor was used to reproduce the original 
critical exponents r) Li and v L ± from the seminal paper [l|]. 

Here the geometric angular factor is determined simply by requiring that only J 2 has the 
same leading singularity as in the standard (p 4 theory. In that case, it is simply given by the 
area of the m-dimensional sphere. 

We can now discuss the fixed points and critical exponents for arbitrary m-axial behavior. 
Note that in the isotropic case, only I2 and J4 have the same leading singularities as in the 
standard (ft 4 field theory. This will lead to a different fixed point for the isotropic behavior, 
as it is going to be shown in a moment. 



VI. THE CRITICAL EXPONENTS FOR THE ANISOTROPIC BEHAVIORS 



A. The dissipative approximation 



The critical exponents were first calculated using this approximation for the uniaxial 
case and the generalization was soon presented for the m-axial case. Details can be found 
in Refs. P0|,P5|]. Here we shall simply quote the results. 

The fixed point at two-loop order is given by: 



6 



u 



S + N 



/4(5iV + 22) 



1 [him ~ 



(2 + N) 
(8 + N) 2 



(174) 



It can be used to obtain the critical exponents i]L2 and ul 2 - 
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+ e 



(2 + TV) 
(8 + N) 2 



'4(5iV + 22) 
(8 + N) 2 




2 + N 
8 + N)' 



(175) 



1 1 2 + N 

UL2 -2 + Z €l 8Tn 



+ l(2 + N) 



+ Ny 



2{UN + 40) [i 2 } m - 2(2 + N) + (8 + N){3 + N) 



(176) 



Using the scaling law j L = ^2(2 — VL2), the exponent Jl is 



1 2 + N 
n = 1 + 2 €L 8TN 

1 (2 + N) 
A (8 + N) 3 



12 + 8iV + N 2 + 4 [i 2 ] m (20 + 7N) 



;i77) 



The evaluation of Feynman diagrams used to obtain these results have some inconve- 
niences as discussed before: the introduction of diagram factors for the integrals I' z and I' b 
(which are divergent when m = 8) is the main trouble. 

The RG analysis by Mergulhao and Carneiro has been used in a attempt to extend the 
calculation for all m. Using the fact that the quartic momenta scale including a and the 
quadratic external momenta scales are equal Diehl and Shpot considered the anisotropic 
problem for general m ]3"T| , |3"2]] . In their first work |HJ, they worked directly in position 
space. After that, using a hybrid approach, going to coordinate or momentum space using 
the free propagator (scaling function) in coordinate space to make the transition according to 
the necessity, they calculated the critical exponents using a minimal subtraction procedure 
which sets the external quartic momenta scales to zero. However, there is a small discrepancy 
among their results for the critical exponents in the cases m = 2, 6 when compared with 
Mergulhao and Carneiro's results using normalization conditions []32]j . For the anisotropic 
cases m ytz d, the exponents are given in terms of integrals to be performed numerically. 
These numerical integrals are meaningful solely if one separates the integration limits on 
the variable v = cr X||X_|_ using the scaling and related functions in the coordinate space 
representation in the integrand up to the maximum value of |v| at \vo\ = 9.3, and replacing 
the asymptotic value of these functions for greater values of v. Note that as the quartic and 
quadratic external momenta are not independent, a cannot be taken dimensionless as done 
by these authors following the invalid argument by Mergulhao and Carneiro. Thus, they 
erroneously concluded that the isotropic case could be encompassed in their expressions for 
the critical exponents in the limit d — » m close to 8. 

Furthermore, this alternative semianalytical method has some drawbacks. First, setting 
the quadratic momenta scale to zero makes impossible the transition from the anisotropic to 
the isotropic case, since the quadratic momenta scale is absent in this case and renormaliza- 
tion group transformations are defined only through the variation of the quartic momenta 
scale. Second, unfortunately the expression of the critical exponents for the anisotropic case 
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are rather cumbersome, given in terms of integrals to be performed numerically. Clearly the 
most convenient answer should present analytical coefficients for each order in e^, in analogy 
to the usual standard (p 4 theory describing the Ising model. 

The calculation of the critical exponents using the dissipative approximation presented 
here has been criticized by Diehl and Shpot |33[] because of the constraint introduced in 
the quartic loop momenta in higher loop Feynman integrals. Despite of this criticism, this 
approximation is in good agreement with recent high-precision numerical data based in 



Monte Carlo simulations for the ANNNI model [34. 22 



B. The orthogonal approximation 

We turn now our attention to the most general approximation for calculating all the 
critical exponents. In order to check the results, we must calculate the critical exponents 
using two different renormalization schemes, namely the normalization conditions and min- 
imal subtraction of dimensional poles. It is going to be shown now that the critical indices 
are the same irrespective of the use of either renormalization prescription. 



1. Normalization conditions and critical exponents 

We defined the bare coupling constants and renormalization functions as 

u OT = u T (l + a\ T u T + a 2T u T ), (178a) 
Z m = 1 + b 2r u\ + b 3r u 3 T , (178b) 

Z^2( T ) = 1 + Ci T U T + C 2t uI, (178c) 

where the constants a iT ,bi T ,c iT depend on Feynman integrals calculated at the convenient 
symmetry points. Depending on the symmetry point, we can calculate critical exponents 
corresponding to correlations perpendicular or parallel to the competing m-dimensional 
subspace. 

The beta-functions and renormalization constants can be rewritten in terms of the con- 
stants defined above in the following form: 

(3 T = -Te L u T [l - a lT u T + 2(a\ T - a 2T )u 2 T ], (179a) 
7<6(r) = -re L u T [2b 2r u T + (36 3r - 2b 2r ai T )u% (179b) 
l^{r) = re L u T [ci T + (2c 2t - c\ T - a lT c lT )u T ]. (179c) 

It is easy to obtain the coefficients above as functions of the integrals calculated at the 
symmetry points. They are given by 

N + 8 r r . , 
air = — a 1 + (180a) 

jy + 8 r (iV + 8) 2 r . 1 (3JV+14) 1 
a 2T = (— ) + [^— W\ m - 24 (18 ° b) 

6 2T = -^^[l + (2[i 2 ] m + ^], (180c) 
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= _ (N + 2) { N + 8) { N + 2)(iV + 8) 1 j_ 
3 1296e| 108e L V 4 L J 48 ; V ; 

Clr = ( ^ +2) [l + Nm6L], (180e) 



(iV + 2)(iV + 5) (^ + 2)^ + 5) ,., 1 
= m + -^—[— r - N m - -]. (180f) 



This is enough to obtain the fixed points at 0(e 2 L ). They are defined by (3 T {u* T ) = 0. 
All the integrals calculated at the symmetry points SP± and SP2 look the same. The 
factor of r = 1,2 drops out at the fixed points, implying that the renormalization group 
transformations performed either over K\ or k 2 will flow to the same fixed point given by 
(u? = u% = u*) 



, (9iV + 42) 



(181) 



The surprising feature of this approximation is that the critical exponents do not depend 
on [ia]mj as the ones obtained using the dissipative approximation. The functions 7^1) and 
7<£2(i) can be written as 

70(i) = ^-^t 1 + ( 2 Nm + ^lH - * ^ '-ul (182) 

(iV + 2) r . . 1 . . . 

70 2 (i) = q u A l + V2\ m tL - (183) 

Replacing the value of the fixed point inside these equations, using the relation among these 
functions and the critical exponents r]L2 and Vl2, we find: 

1 2 N -\- 1 , .6(37V+14) 1 

^ 2 = 2 ei (a7T^F [1 + ^TFTsT " (184) 

1 (iV + 2) l (iV + 2)(iV 2 + 23iV + 60) 2 
" L2 = 2 + 4(A^) eL + 8 (ivT^ ^ (185) 

Notice that the coefficient of each power of is the same as in the pure 4 describing the 
Ising-like behavior. The reduction to the m = case is even simpler using this approxima- 
tion than the reduction to the m = case using the dissipative approximation. Since the 
functions 7^(2) = 27^(1) and 7^ 2 (2) = 27^2 m as a consequence of = 2/3i, we immediately 
conclude that 

, (N + 2) r ,6(3iV + 14) L. . , 

1 , (N + 2) 1 (iV + 2)(iV 2 + 23iV + 60) 2 

4 8(ivT8) eL + 16 (N + W ' (187) 

Thus, at 0(e|), the relation r/ £4 = 2r] L2 is valid. At 0(ej j ), the relation z^ 4 = \vl2 is 
fulfilled. Thus the strong anisotropic scale invariance |M] is exact to the perturbative order 
considered here. The other exponents can be read from the scaling relations. As discussed 
before, they are [o-L2 = «L4 = «l, etc.): 



38 



(jV + 2) (iV + 2)(iV 2 + 22iV + 52) 

7L = 1 + 2(^) eL + 4(]vT8p ^ (188) 

(4-iV) (iV + 2)(iV 2 + 30iV + 56) 2 

ai = w^) eL ww 3 — (189) 

z?- 1 3 (iV + 2)(2iV + l) 

2 2(iVT8) ei + 2(iV + 8)3 ^' (190) 

(A^ 2 + 147V + 60) „ , , 

Note that all these exponents reduce to the Ising-like case when m — 0. In order 
to check the correctness of these exponents, it is convenient to calculate them in another 
renormalization procedure, as we shall see next. 



2. Minimal subtraction and critical exponents 

Usually, in the minimal subtraction renormalization scale, one can have more than one 
coupling, but just one momenta scale, called in most textbooks /i [[36] and named k here. The 
dimensional redefinition performed for the quartic external momenta, allows the description 
of the anisotropic case with two independent momenta scales. The coupling constant has 
two independent flows, induced by K\ and K2. 

If we want to calculate the critical exponents along the competition axes, we set the 
quadratic external momenta perpendicular to the competing subspace to zero. Thus, one can 
introduce the quartic momenta scale Ki in order to compute the normalization functions for 
arbitrary quartic external momenta and demanding that the dimensional poles (logarithmic 
divergences in the momenta) be minimally subtracted. On the other hand, the calculation 
of critical exponents perpendicular to the competing axes can be performed by setting the 
quartic external momenta to zero, introducing K\, calculating the normalization functions 
for arbitrary quadratic external momenta and requiring minimal subtraction. 

In this section, we are not going to calculate explicitly the critical exponents. Instead, 
we are going to calculate the fixed point as well as the functions J<j>( T ) and 7^2 ( T ) at the 
fixed point. As these functions at the fixed point are universal, they should be equal to 
the ones obtained using normalization conditions, leading to the same exponents in either 
renormalization scheme. 

The dimensionless bare couplings and the renormalization functions are defined in min- 
imal subtraction by 

00 

u 0t = u T [l + Ctir(eL)<], (192a) 

i=l 
oo 

%r) = l + EMeL)<, (192b) 

i=l 
oo 

= 1 + £ c -(qX- (192c) 

i=l 

The renormalized vertices 
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r fl(r)(^, Ur, Kr) = Z^T^kr, W r, «t)» (193a) 
r R(r)(^r, Mr, «t) = Zj (r) r [J] (fc ir , U 0t , K t ), (193b) 
r i?(r)(^'^2r,Pr;Mr,/«r) = Z ^2 {t) T^ (k lr , k 2r , p T , U 0t , K T ) , (193c) 

are finite when — > 0, order by order in u T . Note that the external momenta into the bare 
vertices are mutiplied by k~ 1 . Recall that kn = Pi are the external momenta perpendicular to 
the competing axes, whereas ka, = k\ are the external momenta parallel to the m-dimensional 
subspace. The coefficients a ir {ei), b ir {ei) and c iT {ei) are obtained by requiring that the poles 
in tL be minimally subtracted. The bare vertices can now be expressed as 

r$(Av, u 0T , ^ = k 2 T T {l - B 2t u 2 0t + B 3t u 3 0t ), (194a) 

r§(k iT , U 0T , Kr) = < e Mor[l - A 2t U 0t + (A$ + (194b) 
rjj^fclr, k 2 r,Pr] U 0t , «v) = 1 - C lr «0r + (Cff + CgVor" (194c) 

Notice that £?2t is proportional to the integral I3 and _E>3 T is proportional to I5. Note that if 
t = 1, all the integrals should be replaced by their values at zero quartic external momenta. 
In case r = 2, those integrals are calculated at zero quadratic external momenta. 
Explicitly, the coefficients are given by the following integrals: 

A (N + 8) k lT + k 2T T{ k lT + k 3r T{ k 2r + k 3r 

Air = — 7^— [M ) + H ) + H )], (195a) 

A% = {N2 + ™ + 20) jll( kl ^ k2T ) + i|(^±^) + / 2 2 (\ + H (195b) 

4r = ^ [ /4 (~) + 5 Permutations], (195c) 

C lT = ^[J 2 (^±^I) + / 2 (^±^L) + + (195f) 

C$ = + J|(^±^I) + /|(^1±^)], (195g) 

C$ = — — ^ — [h(—) + 5 permutations]. (195h) 

This is sufficient to determine the normalization constants to the loop order desired. 
Requiring minimal subtraction of dimensional poles for the renormalized vertex parts quoted 
above, all the logarithmic integrals in the external momenta appearing in I 2 ,I 3 ,l4, and 1$ 
cancell each other. The result is that the normalization functions and coupling constants 
can be expressed in the form: 



(N + 8) AN + 8) 2 (37V + 14V „ 
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J <t> 2 {r) 



iV + 2 2 (iV + 2)(iV + 8) , (JV + 2)(JV + 8) 1 _g 



144e L 
7V + 2 



1296e| 



5184e L 



1 + ^„ T + [ g±jH^±j)_(^±j) lMr 



6e L ' 36e£ 24e l 

From the renormalization functions one can obtain: 



7<«t) 



7^2 (t) = r- 



(iV + 2) 2 (iV + 2)(iV + 8) 

72 " r 1728 

(iV + 2) 



6 



(196b) 
(196c) 

(197) 
(198) 



The fixed points are defined by f3 T {u* T ) = 0. Then, it is found that the fixed points 
generated by renormalization group transformations over either K\ or are the same and 
is given by: 



6 



u' 



8 + N 



e L {l + e L 



(9N + 42) 
(8 + N) 2 



(199) 



Substitution of this result into the renormalization constants will give at the fixed point 
7<£(t) = Vtj where r\ T are given by Eqs. (184) and (186). In addition, we have 



7>(r) 



(N + 2) 6(iV + 3) 

r— — —e L [l + — T : ^ n e L \. 



(N + 



(N + 



(200) 



This leads to the same exponents v T given in Eqs. (185) and (187), obtained there via 
normalization conditions. Therefore, we have proven the consistency of this picture for 
the anisotropic Lifshitz critical behavior, since the critical indices are independent of the 
renormalization procedure. 



3. Discussion 

The exponent 77^2 obtained here agrees with the calculation performed independently by 
Mukamel P7|. Nevertheless, the exponent r]L4 presented here is at variance with Mukamel's 



37J and, therefore, with the result obtained by Hornreich and Bruce |38| since both works 
agree with each other. 

We are now in position to compare our results with those obtained for the ANNNI 
model in three-dimensional space (e^ = 1.5) representing the uniaxial (m = 1) case using 
Monte Carlo simulations From the numerical viewpoint, there is no sensitive difference 



among the results presented either using the dissipative approximation or the orthogonal 
approximation for the critical exponents perpendicular to the competition axes. Within the 
two significative algarisms precision the exponents using either approximation are given by 
rj L2 = 0.04 and v L2 = 0.73. 

The deviations start in the calculation of 7^2 = 7l- In the dissipative approximation the 
e^-expansion yields 7^ = 1.45. A numerical interpretation has been proposed recently in 
order to improve the results obtained via the e^-expansion when the perturbative parameter 



€l is greater than 1 p2[ . There it was argued that the neglected 0(e L ) could be relevant to 
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the calculation of, say, 7^. The basic idea is to replace the numerical values of vl2 and r]L2 
directly into the scaling laws in order to obtain the other critical exponents. In this way one 
obtains 7 L = 1.43, a L = 0.18 and (3 L = 0.20. 

On the other hand, using the e^-expansion results for j L , and (3l obtained via the 
orthogonal approximation one finds jl = 1-42, a L = 0.05 and (3 L = 0.26. The numerical 
"ansatz" described above gives again jl = 1.43, ol = 0.18 and fit, = 0.20, since t}l 2 and v^ 2 
have the same numerical values in either approximation. These numbers should be compared 
with the newest Monte Carlo simulations output, namely, 7^ = 1.36±0.03, = 0.18±0.02 
and (3 L = 0.238 ±0.005. 

Thus, the greater the mean field values for the exponents, the better are the their numer- 
ical values using the e^-expansion when is not a small number. Otherwise, the numerical 
"ansatz" yields a rather good agreement with the Monte Carlo results, as in the case for the 
exponents and /?£,. This shows that the new results displayed here are consistent with 
the best numerical values available for the ANNNI model. 



VII. THE CRITICAL EXPONENTS FOR THE ISOTROPIC SYSTEMS 

As the isotropic behavior presents just one external momenta scale, its analysis is simpler 
than the one used to describe the anisotropic behavior, where two external momenta scales 
are present. Besides, the only manner to attack this problem is to use the orthogonal 
approximation, for the dissipative approximation does not work as it was discussed before. 



A. Critical exponents in normalization conditions 

The bare coupling constants and renormalization functions are defined as 

^03 = u 3 (l + a 13 u 3 + a 23 ul), (201a) 
Z m = 1 + b 23 u 2 3 + b 33 ul (201b) 
Z^2 {3) = 1 + c 13 u 3 + c 23 ul, (201c) 

where the constants 0*3,613, 0*3 depend on Feynman integrals calculated at the symmetry 
point named henceafter SP 3 . Only the external momenta scale k 3 parallel to the competing 
m-dimensional subspace arises in this isotropic case. 

The beta-function and renormalization constants are written in terms of the constants 
defined above in the following manner: 

P 3 = -e L u 3 [l - a 13 u 3 + 2(a 2 13 - 023)^3], (202a) 
7^(3) = -e L u 3 [2b 23 u 3 + (36 33 - 2623013)^3], (202b) 
1^(3) = e L u 3 [c 13 + (2c 23 - Ci 3 - a 13 c 13 )u 3 ]. (202c) 

The coefficients above are obtained as functions of the integrals calculated at the sym- 
metry point. They read 
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013 



023 



>23 



'33 



Cl3 



C23 



,iV + 8, 9 r 2iV 2 + 23N + 86, 



6e L 

(iV + 2 



+ [ 
[1 + 



144e, 



(iV + 2)(iV + 8) 
20736e L 



288e L 

(iV + 2)(^ + 8 
25924 

(N + 2)(N + h) (N + 2)(2N + 7) 



36 e | 



The fixed point is denned by 03(uQ 

6 



Mo 



8 + N 



CL 1 



+ 144e L 
0. Therefore, it is given by 
1 



(203a) 
(203b) 
(203c) 
(203d) 
(203e) 
(203f) 



6L 2 



1 (9N + A2) 

2 + (8 + N) 2 



(204) 



Note that this fixed point is different from that appearing in the anisotropic behavior and 
cannot be obtained from it in a smooth way The functions 7^(3) and 7</> 2 (3) can be written 

as 



7^(3) = 

7^(3) 



{N + 2) H 5 . 

■[1+ o 6 ^ 



144 
(iV + 2) 

6 



(N + 2)(N + 



3456 



-u 
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n 1 1 1 

%[1 + -e L - -%]• 



(205) 
(206) 



Replacing the value of the fixed point inside these equations, using the relation among these 
functions and the critical exponents 77^4 and ula, we find: 



VL4 ~ 4 €l JN + 



N + 2 u ,3(3iV + 14) 1. 



(N + 



VhA 4 + 16(N + 8)' 



(N + 2) 1 (Ar + 2)(N 2 + 23^ + 60) 2 



256 



(N + 



(207) 
(208) 



These exponents are different from those originally obtained in Ref. |L|. The coefficient of the 
e| term in the exponent 77^4 is positive, consistent with its counterpart in the anisotropic 
cases as well as in the Ising-like case. One learns that only the quartic momenta is not 
sufficient to induce its change of sign. The exponent v Li agrees at 0(e/J with that presented 
in Ref. |I| but naturally disagrees at 0(e 2 L ), since it depends on the value of 771,4 at 0(e 2 L ). 
Besides, the critical index 77^4 is obtained at 0(e 3 L ) here for the first time. 

Now using the scaling relations derived for the isotropic case we obtain immediately 



7L4 = 1 + 



(N + 2) 



4(iV + 
(4 -AO 



e L + 



(N + 2)(N 2 + 19N + 28) 



4(N + 8] 



64(AT + 8) 3 
(N + 2)(N 2 + 9N + 68) 
32(N + 8) 3 



(209) 
(210) 
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1 3 (iV + 2)(iV 2 + /V + 108) ^ 

^ L4 "2 4(iV + 8) L 64(iV + 8)3 ^ 

1 (jy2 + 147V + 60) 2 
^ 4 = 3+ 2 £i+ 8(iV + 8) 2 ^ (212) 

These exponents are obtained here for the first time at 0(e 2 L ). Formerly the lack of a set 
of scaling laws for the isotropic case did not allow these findings. In order to check these 
results, let us analyse the situation using the minimal subtraction scheme. 



B. Critical exponents in minimal subtraction 

We proceed analogously as in the isotropic case. We just replace the subscript r = 3 and 
keep in mind that the Feynman integrals are calculated in the isotropic case d = m close to 
8. Minimal subtraction of dimensional poles in the renormalized vertex T^ 4 ^ implies that 
the bare dimensionless coupling constant can be expressed in the form: 



M (iV + 8) AN + 8) 2 (37V + 14), ,, 



The fixed point can be easily found to be 



6 9(3iV+14) 
U * = [N^) L + (N + 8)3 ^ (2M) 



The normalization constants are given by: 

, (# + 2)..2 



(N + 2KAT + 8) (7V + 2)(7V + 8)- 
2592e| + 20736e L 



~ (N + 2) 

Z<f> 2 (3) — H ~ M 3 

6e L 

The functions 7^(3) and 7^(3) are given by the following expressions: 

(N + 2) , (iV + 2)(/V + 8) 3 
70(3) = ggig «s, (217a) 

(iV + 2), 1 9x , 
7^(s) = k -g-^(«s- i «g). (217b) 

Using these results the function 7^3) at the fixed point yields the value of 77^4 as obtained 
in (207), whereas the function 7^2(3) at the fixed point reads 

(7V + 2) 3(iV + 3) 

7 ^ 3 ) = (iV^ eL[1 + 7wW i] ' ( } 



44 



which is the same as that obtained in the fixed point using normalization conditions and leads 
to the same critical exponent zv L4 from (208) as the reader is invited to check. Therefore, 
the complete equivalence between the two renormalization schemes is assured. 

Notice that the critical exponent rj^ for the isotropic case is different from the original 
result Ref. J1J. Since we have checked our results using two distinct renormalization schemes 
as shown above, the critical indices presented by those authors should be checked using more 
than one renormalization procedure in order to clarify this discrepancy. 



VIII. CONCLUSIONS AND PERSPECTIVES 

All the critical exponents for the m-axial Lifshitz critical behavior for the anisotropic 
(1 < m < d — 1) and the isotropic (d = m close to 8) cases are explicitly derived at 0(e 2 L ). 
We have shown that up to the loop order considered in this work strong anisotropic scaling 
theory holds since the relations v L ± = \vl2 and rj Li = 2r] L2 are exact. The exponents 
associated to critical correlations perpendicular to the competing axes easily reduces to the 
Ising-like exponents when m = 0, the only difference being the perturbation parameter ei 
replacing the usual e in noncompeting systems. These relations imply that the new scaling 
laws, obtained here using two independent renormalization group transformations, reduce 
to the ones previously found by Hornreich, Luban and Shtrikman [|TJ. 

Moreover, to our knowledge this is the first time that all the exponents for the isotropic 
behavior are obtained explicitly through the use of the new scaling relations presented |52" 



Besides, they are shown explicitly not to be recoverable from the anisotropic situation in the 
limit d — > m. The structure of the Feynman integrals in the isotropic case indicates that it 
deserves a especial treatment when compared with the anisotropic situation as clarified in 
this article. 

The new results for the calculation of arbitrary loop Feynman integrals are obtained 
by demanding that they are homogeneous functions of arbitrary external momenta. Even 
though the calculations are carried out in a given order in perturbation theory, the author 
is convinced, however, that the conclusions hold to all orders. 

The simple analytical expressions for each coefficient in the e^-expansion of the critical 
indices are rather encouraging to proceed the evaluation of other universal amounts, like 
critical amplitudes |55| . It would be interesting to compare some experimental results avail- 
able for MnP like the specific heat critical amplitude ratio [|l^] with theoretical calculations 
within the context of an e^-expansion using the techniques described in the present work. 
In addition, a thorough RG analysis to prove that all amplitude ratios are indeed universal 
for the Lifshitz critical behavior was not done yet. Actually, the idea presented in this work 
might be suitable to demonstrate the universal character of the above mentioned critical 
ratios and calculate all of them. 

Other problems can be pursued using the present method. The treatment of finite-size 



effects for the Lifshitz behavior can be devised in analogy to the noncompeting situation |39], 
|4"0"f . The systems may be finite (or semi-infinite) along one (or several) of their dimensions, 
but they are of infinite extent in the remaining directions. Examples include systems which 
are finite in all directions, such as a (hyper) cube of size L, and systems which are of 
infinite size in d' — d — 1 dimensions but are either of finite thickness L along the remaining 
direction (rf-dimensional layered geometry) or of a semi-infinite extension. The presence of 
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geometrical restrictions on the domain of systems also requires the introduction of boundary 
conditions (periodic, antiperiodic, Dirichlet and Neumann) satisfied by the order parameter 
on the surfaces. In particular, the validity limits of the e^-expansion for these systems and 



the approach to bulk criticality in a layered geometry can be studied [41 



Recently, typical surface phenomena in noncompeting systems were generalized to corn- 
However, as far as 



peting systems using Monte Carlo simulations for the ANNNI model [42 



the Lifshitz behavior is concerned, a theoretical description of these systems is still lacking. 
The field-theoretical framework just presented might be useful to address this problem. 

The quest towards a generalization of the Lifshitz universality classes whenever arbitrary 
momenta powers arise in the Lagrangian (1) as the effect of further competition is quite a 
fascinating issue |43| . It is expected that it can be solved along the same lines described in 
this work [[44]]. 

In summary, we have described the Lifshitz critical behavior in its complete generality 
in what concerns its critical exponents. We have presented new field theory renormalization 
group methods which resulted in new analytical expressions for all the critical indices in the 
isotropic as well as in the anisotropic cases at least at 0(e 2 L ). We hope our findings will be 
useful to unveil further issues related to the physics of competing systems. 
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